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SUMMARY 


This  final  report  on  Contract  F30602 -80-C-0042  summar¬ 
izes  the  AIL  investigation  to  determine  and  demonstrate  suitable 
algorithms  and  correspondi ng  computer  software  in  order  to  pro¬ 
vide  open-loop  adaptive  control  of  the  Microwave  Transversal 
Equalizer  (MTE)  previously  developed  by  AIL  under  Contract  F3062- 
78-C-0352. 

The  program  task  objectives  have  been  successfully  accomp¬ 
lished  by  AIL  with  the  application  of  Fast  Fourier  Transform  ( F FT ) 
techniques  to  provide  an  algorithm  programmed  on  the  RADC  HP  2100A 
computer.  This  procedure  will  ultimately  be  utilized  to  determine 
the  MTE  amplitude  and  time  delay  adjustments  necessary  to  equalize 
the  distortion  introduced  by  an  arbitrary  network  in  series  with 
the  equalizer. 

Verification  of  the  developed  FFT  and  associated  software 
program  was  accomplished  at  AIL  with  a  DEC-20  computer  using 
artificially  generated  simulated  time  sidelobe  distortion  typical 
of  the  expected  range  of  levels.  Translation  of  the  algorithm 
from  the  DEC-20  computer  format  to  the  HP  2100A  computer  format 
has  been  completed,  and  verification  of  the  software  program  using 
RADC  I  and  Q  data  is  underway.  Preliminary  data  suggests  that  the 
program  translation  task  has  been  successfully  accomplished  by  AIL. 
It  is  expected  that  optimal  data  acquisition  techniques  and/or 
preferred  computer  interface  procedures  will  be  defined  as  a 
result  of  related  pos t -del i very  activities  at  RADC. 
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1.  INTRODUCTION 


The  objective  of  Contract  F30602-80-C-0042  is  to  develop 
an  algorithm,  and  relevant  computer  software  programming  for  the 
RADC  HP  2100A  computer,  to  facilitate  open-loop  control  of  the 
MTE  in  order  to  minimize  time  si  delobe  distortion.  The  result¬ 
ant  algorithm  using  FFT  techniques  and  associated  programming 
has  been  validated  with  artificially  generated  simulated  distor¬ 
tion  data  using  the  AIL  DEC-20  computer  and  the  HP  2100A  computer 
at  RADC.  A  similar  post-delivery  validation  effort  will  be  con¬ 
ducted  by  RADC  personnel  using  real-time  distortion  data  and  tne 
HP  2100A  computer. 

Section  2.0  will  contain  a  brief  discussion  of  MTE  opera¬ 
tional  fundamentals  in  order  to  establish  an  inter-relationship 
with  the  algorithm  development  program.  Section  3.0  will  pre¬ 
sent  a  detailed  treatment  of  the  FFT  techniques  applied  to 
algorithm  development  for  the  MTE.  Section  4.0  will  discuss 
the  resultant  computer  program  including  applicable  operating 
and  test  procedures.  Finally,  Section  5.0  will  present  con¬ 
clusions  and  recommendations  for  a  logical  continuation  of 
effort  required  to  provide  closed-loop  operation  of  the  MTE. 
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2.  MTE  TECHNICAL  DISCUSSION 


A  discussion  of  distortion  considerations  and  the  funda 
mentals  of  the  MTE  operation  are  presented  in  order  to  provide 
the  technical  basis  for  computer  control  and  adaptive  equaliza¬ 
tion  of  resultant  time  sidelobe  distortion.  A  detailed  descrip¬ 
tion  of  the  MTE,  including  operating  procedures,  is  included  in 
Final  Report  RADC-TR-80-1 21  of  31  January  1980  titled  "Microwave 
Transversal  Equalizer". 


2.1  PAIRED  ECHO  THEORY  AND  DISTORTION  CONSIDERATIONS 

Algorithms  for  determining  the  settings  of  attenuation 
and  phase  for  each  loop  of  a  transversal  equalizer  rely  heavily 
on  the  paired  echo  concept.  Accordingly,  a  brief  review  of  this 
theory  is  given  below  stressing  its  applicability  to  the  present 
program . 


2.1.1  Paired  Echo  Theory 

The  paired  echo  concept  and  its  application  to  the 
design  of  a  transveral  equalizer  (MTE)  was  described  in  Ref.  1 
It  should  be  noted  that  this  work  w^s  supported  by  RADC.  This 
paper  also  gives  a  bibliography  of  prior  work  on  non-microwave 
MTE's  and  the  now  classic  reference  on  Chirp  radar  by  Klander, 
et  al .  The  latter  emphasized  the  use  of  paired  echo  theory  to 
predict  time  sidelobe  levels  in  Chirp  radars.  The  essential 
aspects  of  the  paired  echo  theory  are  summarized  below. 


Ref.  1  -  J.J.  Taub  and  G.P.  Kurpis,  "Microwave  Transversal  Equal 
izer",  Microwave  Journal,  1969. 
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Signals  which  can  be  characterized  by  a  time  function 
of  its  Fourier  transform  (amplitude  and  phase  vs  frequency  func¬ 
tions),  are  subject  to  distortion  when  propagating  through  a  micro- 
wave  component  or  a  system  of  microwave  components  (such  as  a  Chirp 
radar).  This  distortion  occurs  because  the  system's  transfer  func¬ 
tion  possesses  neither  perfectly  constant  gain  nor  perfectly  linear 
phase  over  the  spectrum  of  the  signal. 

We  typically  represent  the  transfer  function  in  the  fre¬ 
quency  domain  and  need  to  predict  its  effect  on  time  domain  dis¬ 
tortion.  By  using  the  paired  theory  we  can  make  relatively  quick 
conversions  between  the  frequency  domain  and  the  time  domain  and 
vice  versa. 

The  system  transfer  function,  or  frequency  response,  of 
an  aribtrary  system  can  be  defined  as: 

H(u>)  *  A(u)  ejB^  (2-1) 

where  A  and  B  are  respective  gain  and  phase  functions  defined  in 
Reference  1 . 

We  must  now  assume  that  the  signal  transmitted  through 
the  system  has  spectral  components  that  are  band  limited.  This 
is  a  safe  assumption  for  most  systems.  For  example,  in  many  radar 
applications  the  signal's  spectrum  is  confined  to  the  electronic 
bandwidth  of  the  final  power  amplifier.  Within  a  band  limit  of 
to  we  can  rigorously  represent  the  A  and  B  functions  as 

Fourier  series: 
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In  a  distortionless  network  A(^)  =  aQ  and  B(^)  =  b0w' 

Hence  the  remaining  terms  contribute  to  distortion.  This 
representation  leads  to  the  time  doma in  response  including 
distortion  and  so  that  for  a  band  limited  signal 
V.(t-)  and  its  time  function  equivalent  V ^  ( t )  we  obtain 

V *>  ■  ao  #t  <T  +  b0>  *  *„  £,  E„*Vi  <t+l>o  -  "*> 

n  =  l 

N 

*  ao  „E,  E„-  V,  +"') 

where  V  ( t )  is  the  output  time  function  and  the  constants  E  .  and 
o  n  + 

Ep_  represent  echo  or  time  sidelobe  amplitude  levels.  These  l^pl  s 

are  related  to  the  Fourier  distortion  coefficients  in  equations 

(2-2)  and  (2-3).  For  the  case  where  0  and  =  0  they  reduce  to 

o  o 
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(2-5) 


A  typical  example  of  a  distorted  microwave  pulse  is  given 
in  Figure  2-1.  The  echoes  are  clearly  displayed.  This  analysis 
desmonstrates  that  a  frequency  domain  transfer  function  can  be 
represented  by  a  Fourier  series  which  enables  rapid  calculation 
of  the  time  domain  response.  Conversely,  a  measurement  of  the 
time  domain  response  which  yields  si  delobe  levels  can  be  used  to 
rapidly  calculate  the  a  and  b  coefficients  thereby  yielding  the 
frequency  domain  transfer  function.  The  MTE  settings  cancel 
distortion  by  injecting  equal  and  opposite  echoes  in  cascade 
with  the  system. 

2.1.2  Use  of  Paired  Echo  Theory 

Since  the  MTE  design  and  adjustment  procedures  are  based 
on  cancelling  distortion  echoes,  algorithms  can  be  developed  by 
taking  either  frequency  domain  (I  and  Q  data  )  measurements  or 
time  domain  responses  (time  sidelobe  levels  in  d3)  and  convert¬ 
ing  them  into  the  necessary  MTE  loop  attenuation  and  phase  set¬ 
tings  to  produce  cancelling  echoes.  The  key  point  of  this  dis¬ 
cussion  is  that  use  of  paired  echo  theory  simplifies  computation 
and  therefore  significantly  reduces  computer  time.  Furthermore 
it  provides  the  flexibility  to  work  with  either  frequency  or 
time  domain  data.  In  this  program  the  actual  procedure  used 
was  to  use  I  and  Q  data  (frequency  domain  measurements). 
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MAINLOBt 


Figure  21.  Time  Sidelobe  Distortion. 


2.2  MICROWAVE  TRANSVERAL  EQUALIZER  OPERATION 

The  MTE  has  been  designed  by  AIL  to  cancel  the  para¬ 
sitic  time  sidelobes  caused  by  distortion,  by  generating  its 
own  artificial  equalizer  sidelobes,  whose  envelopes  are  prop¬ 
erly  delayed  (or  advanced)  to  coincide  in  time  with  the  para¬ 
sitic  sidelobes.  In  addition,  the  amplitude  and  phase  of  each 
artificial  sidelobe  can  be  separately  adjusted  to  produce  an 
equal  amplitude  phase  reversed  replica  of  the  corresponding 
parasitic  side! obe . 

Generation  of  the  equalizer  echoes  is  accomplished 
by  tapping  off  portions  of  energy  from  the  main  signal  lobe 
on  the  main  transmission  path,  directing  these  portions  through 
properly  adjusted  delay  lines  and  reintroducing  the  delayed 
signal  portions  back  into  the  main  transmission  path  with 
proper  amplitudes  and  polarity.  The  time  domain  output 
response  of  an  ideal  MTE  is  shown  in  Figure  2-2.  It  should 
be  noted  however  that  the  actual  response  will  consist  of  a 
train  of  echoes  with  decreasing  amplitudes  (0.75  dB/tap)  due 
to  the  insertion  loss  of  the  mainline  tapped  delay  elements. 

A  correction  for  the  aforementioned  MTE  echo  transfer  char¬ 
acteristics  will  be  necessary  in  order  to  achieve  the  desired 
tap  attenuator  setting  with  adaptive  closed  loop  operation 
(computer  controlled). 
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2-2.  Time  Domain  Output  of  Microwave  Transversal  Equalizer. 


The  MTE  consists  of  four  cascade  connected  microwave 


integrated  circuit  mainlines  and  32  secondary  lines  for  the 
purpose  of  generating  a  train  of  8  pre-  and  post-mai nl obe  echoes 
(16  total).  A  3.3  nsec  time  delay  between  echoes  of  each  train 
is  provided.*  The  amplitude  level  and  time  delay  of  the  selected 
16  of  the  32  possible  echoes  are  individually  adjusted  with  the 
PIN  diode  variable  attenuator  (electronically  controlled),  a 
line  stretcher  variable  delay  (mechanically  controlled),  included 
in  each  secondary  line.  The  remaining  16  echo  taps  are  terminated, 
however  these  taps  are  available  for  use  as  required. 

It  should  be  noted  that  the  transfer  characteristics  of 
each  of  the  16  variable  attenuators,  as  well  as  the  16  variable 
solid  state  time  delay  units  (when  they  are  developed)  are  required 
to  properly  achieve  adaptive  MTE  control. 


3.  ALGORITHM  DEVELOPMENT 

An  algorithm  has  been  developed  using  FFT  techniques 
to  provide  open  loop  control  of  the  MTE.  The  technique  will 
be  described  in  terms  of  the  relevant  input  and  output  formats, 
as  well  as  the  FFT  algorithm. 

3.1  INPUT  FORMAT 

The  input  to  the  FFT  consists  of  a*set  of  complex  correl¬ 
ations  as  a  function  of  frequency.  These  samples  can  be  expressed 
by 

H( jcu1 )  *  I  (w,.)  +  j Q ( u)  1  )  (3-1) 

where  . 

w.  is  the  ic  frequency 

I(u)1 )  is  the  in-phase  frequency  response 

Q(to^)  is  the  quad-phase  frequency  response 

For  the  problem  at  hand,  the  following  frequencies  are  considered: 
U>1  =  27 T  [fL  +  ( i  - 1  )  Af],  i  =  1  ,2  , .  256  (  3-2) 

where 

fL  =  3.1  GHz 

Af  =  1  .1  7647MHz  -  1.18  MHz 

The  last  frequency  is  found  by  setting  i  =  256;  that  is, 
fH  *  3-1  +  255  x  1.18  x  10'3  =  3.4  GHz 

As  a  result,  a  total  of  256  pairs  of  measurement  data  are  involved 
for  a  bandwidth  of  300  MHz  as  shown  in  the  format  presented  in 
Tabl e  3-1 . 
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TABLE  3-1.  TAPE  FORMAT  FOR  THE  FREQUENCY  RESPONSE 


Word  Index 

Notation 

Frequency 

1 

Kj*!  ) 

3.1  GHz 

2 

r 

Q(  jw1 ) 

3.1  GHz 

3 

I ( jw  2 ) 

3.1  +  Af 

4 

C\J 

3 

O' 

3.1  +  Af 

0 

0 

0 

0 

0 

0 

0 

0 

0 

509 

I(j»255) 

3.1  +254  At 

510 

Q(jw255) 

3.1  +  254  Af 

51  1 

I(J“  2s<! 

3.4  GHz 

512 

25 i 

3.4  GHz 

Note : 

1  . 


-  3.4  -  3.1 
’  255  x  io’3 


1.18  MHz 


2.  Each  word  is  represented  by  8  bits  in  binary  2's  complemert 
notation  ;  i . e . : 

-128  <  I,g  <  +  127 


The  above  representation  defines  the  formula  transforma¬ 
tion  required  from  the  input  medium  to  the  computer. 
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The  input  samples  are  represented  by 
an  8  bit  binary  two's  complement  notation.  Thus,  a  word  is 
defined  to  be  8  bit  with  the  value  of  the  word  ranged  from  127 
to  -128.  In  floating  point  notation  the  value  is  ranged  from 

127/128  =  0.992188  to 
-128/128  =  -1 

3.2  OUTPUT  FORMAT 

The  output  of  the  FFT  consists  of  a  set  of  =  17 
complex  MTE  tap  coefficients  which  represents  a  truncated  Fourier 
series  representation  of  the  equalizer  transfer  function.  The 
total  number  of  relative  equalizer  taps  is  sixteen  with  the  center 
tap  (number  9)  taken  as  the  reference  tap. 

The  output  tap  coefficient  will  be  presented  by  decimal 
representation  with  the  absolute  magnitude  of  the  tap  less  than 
unity,  that  is:  -1  <  <  1  for  i  =  1,2, ...17  and  phase  angle  is 

expressed  in  degree  units. 

Further  data  format  transformation  will  be  needed  to 
match  the  exact  setting  on  the  MTE  attenuator  dial. 

3.3  FFT  ALGORITHM 

Consider  again  the  measured  frequency  response  of  a  net¬ 
work  with  a  transfer  function  H ( jou ^ ) .  The  measurement  can  be  func¬ 
tionally  described  by  the  block  diagram  shown  in  Figure  3-1.  The 
input  signal  to  the  device  can  be  written  as 

X.(t)  =  A,-  cos  (u).t  +  0o)  (3-3) 
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I(jw)  Q(jw) 


FIGURE  3-1.  MEASUREMENT  SETUP 
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where 


.  is  the  it^1  in-band  frequency 
0Q  is  the  input  phase  angle 
A  is  the  amplitude  of  the  signal 

The  response  of  the  network  is 

Y  .j  ( t )  =  Bi  cos  (oj^t  +  9.j }  (3-4) 

where 

9 ^  is  the  corresponding  phase  angle  for  the  ith  signal 
is  the  corresponding  amplitude  for  the  i**1  signal 

The  quadrature  hybrid  is  used  to  generate  the  in-phase  and  the 
quad-phase  components: 


I  (  ju)  .  ) 

1 

T 

;T  a.b. 

cos 

(^it+0o) 

cos  (coit+9i ) 

dt 

(3-5) 

and 

Q(  Ja>i  ) 

1 

T 

0 

T 

/  A.B. 

si  n 

(»,t+0o) 

cos  (a) .  t+9 ^  ) 

dt 

(3-6) 

o 


where  T  is  defined  as  the  correlation  time, 
can  be  reduced  to 

AB. 

I( jui )  =  -f-  cos  ( 0o-Oi ) 

AB  ■ 

Q ( joa i )  =  ~2~  sin  (0o'9i  * 


The  correlator  output 

(3-7) 

(3-8) 


The  gain  and  phase  response  as  a  function  of  frequency  can  be 
expressed  by 


Gi  ■  /  if  +  if  *  \  •  B,  <3-9> 

.  q . 

9.  =  tan"1  (y1)  +  0  (3-10) 

1  1  i  0 
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where  the  short  hand  notation  is  adapted;  that  is 
I,-  =  I(jwi )  and  Qi  3  0 ( j ai ^ ) 

Since  A  and  are  not  a  function  of  frequency,  these  measurement 
data  represent  the  transfer  function  H ( joj ^ ) .  The  impulse  response 
h(nT)  can  be  found  by  discrete  Fourier  transform  method. 

A  few  assumptions  will  be  made  such  that  the  Fourier  series 
approach  can  be  applied: 

(1)  Signals  passing  through  the  network  will  be  bandlimited  to 
300  MHz.  To  provide  equalization  over  this  bandwidth  requires 

a  corresponding  time  delay.  This  defines  the  delay  between  taps 
to  be: 

T  *  a—1  f  -  3  - - - r  =  3.33  nsec 

TH"TL  300  x  10 

(2)  Since  any  digital  filter  spectrum  to  be  realized  will  be 
made  periodic  of  period  f  ,  then  the  spectrum  can  be 
represented  by  a  Fourier  series: 

oo 

H(jco)  3  E  ct  exp  ( j ku>T )  (3-11) 

k=  —  K 

where  ck  is  the  kth  Fourier  coefficient  defined  by 

ws 

Cn  3  ~ —  /  2  H ( ju» )  exp  (-jkuiT)  du  (3-12) 

*  '“S  „s 

'  1 

Generally,  if  the  gain  response  is  an  even  function  and  the  phase 
response  is  an  odd  function,  that  is 

G(ju>)  3  G(  -  ju> )  (3-13) 

9(J«)  3  -Q(-jw)  (3-14) 
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then  the  resultant  coefficients  will  be  real.  Note  that  these 
requirements  that  G ( j to )  and  9(jw)  be  even  and  odd  respectively, 
lead  to  similar  requirements  that  the  in-phase  and  the  quadrature-phase 
components  to  be  even  and  odd,  respectively.  In  practice,  however, 
these  symmetries  cannot  be  assumed  as  shown  in  Figure  3-2.  As  a 
result,  the  coefficients  will  be  complex;  that  is: 


ck  =  ak  +  j  b  ,  k  »  1,2 . 17  (3-15) 

where 

9k  ■  /  4  *  *1  (3-16) 

9  k  =  tan"1  (bk/ak)  (3-17) 


To  achieve  a  set  of  real  equalizer  coefficients  usually 
means  to  double  the  upper  frequency  bound  of  the  signal  samples  (to  )  . 
When  the  upper  frequency  is  doubled,  the  second  half  of  the  spectrum 
is  copied  from  the  original  bandwidth  with  the  proper  polarity 
attached.  Examples  to  illustrate  these  properties  will  be  presented 
in  the  next  section. 

3.4  FFT  ILLUSTRATIVE  EXAMPLES 

A  few  examples  illustrates  the  Fourier  transform  approach 
to  well  defined  H(jw)  will  be  given  in  this  section. 

3.4.1  Lowpass  Filter 

Figure  3-3  shows  a  specification  for  a  lowpass  transfer 
function  with  cutoff  frequency  set  at  0.5  rad/sec  and  the  upper  fre¬ 
quency  bound  (u>  }  set  at  2  rad/sec.  The  spectrum  between  u^/2  to 
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FIGURE 


(a)  GAIN  SPECTRUM 

e(j<o) 


3-2.  GAIN  AND  PHASE  RESPONSE  (GENERAL  CASE  -  ASYMMETRIC) 
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gjs  is  duplicated  to  provide  an  even  function  for  the  in-phase 
response  I(jw).  The  quad-phase  response  is  set  to  zero. 

The  Fourier  coefficients  can  be  readily  found  to  be 

ck  »  {—  sin  (kj,  k  =  -N,  -N+1....N 

k  f  0  and  N  =  8 


and  CQ  =  0.5. 

For  the  case  N  =  8,  or  17-tap  equalizer  the  response  of 
the  network  can  be  plotted  as  shown  in  Figure  3-4. 

Generally,  the  performance  of  such  a  filter  indicates 
a  passband  ripple  of  0.75  dB  and  a  stopband  rejection  of  approxi¬ 
mately  -20  dB.  Direct  truncation  of  the  Fourier  series  leads  to 
the  well-known  Gibbs  phenomenon.  The  ripple,  in  generaly,  cannot 
be  reduced  by  simply  including  more  taps  as  displayed  by  the  25 
tap  version  shown  in  Figure  3-5. 

The  reduction  of  the  passband  ripple  and  stopband  attenu¬ 
ation  is  later  approached  by  finding  a  time-limited  function  whose 
Fourier  transform  best  approximates  a  bandlimited  function.  This 
approach  leads  to  the  well-known  Kaiser  window  expressed  by: 


I0  (B/l  -  (K/N)2) 

^TbT 


|(K)  = 


■N  <  K  <N 


where  8  is  a  constant  (1  <  8  <  10) 

IQ  is  the  modified  Bessel  function  of  order  zero 
N  is  half  the  number  of  taps 

The  constant  g  can  be  determined  experimentally.  Figure  3-6  shows 
the  same  17  taps  frequency  response  with  Kaisier  window.  The 
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FIGURE  3-4.  LOWPASS  FILTER  RESPONSE  (17  Taps)  DIRECT  TRUNCATION 


FIGURE  3-5.  LOWPASS  FILTER  RESPONSE  (25  Taps)  DIRECTED  TRUNCATION 
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FIGURE  3-6.  17  -  Tap  LOWPASS  FILTER  WITH  WINDOW 


passband  ripple  is  essentially  vanished  and  the  stopband  rejec¬ 
tion  Is  reduced  from  -20  dB  to  -38  dB,  Note  that  the  windowing 
operation  is  a  proper  scale  function  applied  to  the  impulse 
response  of  the  filter.  It  does  not  change  the  hardware  con¬ 
figuration  of  the  filter. 

3.4.2  A  Differentiator 

The  response  of  a  differentiator  can  be  written  by 
the  equation 

H(jui)  *  jw,  for  -oj s / 2  <  to  <  Wj/2 

This  is  the  case  that  the  in-phase  response  is  zero  and  the  quad 
phase  Is  an  odd  function.  As  a  result,  the  coefficient  (or  the 
Impulse  response)  is  real  and  odd  functions;  that  is, 

-  ws 

c.  =  —  /  2  w  sin  (kuiT)  dui 

K  <i)_ 

S  0 

1  27T 

-  f  a)  sin  (kuT)  dm  where  T  *  —  =  tt  ,  w  *2 
o  us  s 

■  -  k  =  -N,  -N+1  ,  -1,1,  ...N 

k  =  0 

where  k  =  0 ,  CQ  =  0 . 

Figure  3-7  Is  the  response  of  the  filter  with  17  taps  and  Figure 
3-8  is  the  response  with  Kaiser  window.  The  parameter  B  for  the 
Kaiser  window  is  set  to  3  to  reduce  the  ripple  due  to  direct 
truncation  of  the  Fourier  series  at  N  *  8. 
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IDEAL 


FIGURE  3-8.  (17  Tap)  DIFFERENTIATOR  RESPONSE  WITH  WINDOW 
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Hilbert  Transformer 


The  frequency  response  of  a  Hilbert  transformer  can 
be  defined  as 


H  ( j  u) ) 


0  <_  co  <  cos  /  2 
10  s 

-^  <  a,  < 


The  discrete  Fourier  transform  can  be  expressed  by  the  impulse 
response 


.  /  \  \  s  i  n  i  mr/  2 ) 

"<">  *  >  ~  (mr/2)  ' 


n  f  0 
n  =  0 


The  comparison  between  a  17  tap  filter  and  the  ideal  transform  is 
shown  in  Figure  3-9. 

3.4.4  Equalizer  for  an  Arbitrary  Function 

The  transfer  function  of  an  arbitrary  frequency  response 

H(j(0)  =  I  (  J03  )  +  Q  ( j(0  ) 

over  the  band  of  interest  is  considered.  Since  the  samples  are 
not  symmetrical  in  any  sense,  the  resulting  tap  coefficients  will 
be  complex  and  asymmetrical.  The  matching  of  the  ideal  versus 
17  tap  equalizer  is  plotted  in  Figure  3-10  and  the  resulting  com¬ 
plex  taps  are  tabulated  in  Table  3-2.  Note  that  the  equalizer  taps 
are  complex: 

ck  =  Xk  +  j  Yk’  k  =  1  ’2 . 17 
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RESPONSE  OF  A  {17  Tap)  HILBERT  TRANSFORMER  (EXAMPLE  3) 


and  the  gain  and  phase  adjustment  can  be  found  by  the  transforma¬ 
tion 

gk  *  10  x  log  (Xk^  +  Yj^) 

0k  *  tan-1  (Yk/*k) 


TABLE  3-2.  RESULTING  EQUALIZER  TAPS 


Index 

Xi 

Y  . 

1 

0.438043 

-0.036516 

2 

0.931036 

0.043327 

3 

-0.092596 

-0.064665 

4 

0.086571 

-0.686610 

5 

-0.074379 

-0.573472 

6 

0.132042 

-0.643479 

7 

-0.213581 

-0.047679 

8 

3.133260 

0.643330 

9 

2.183238 

4.936495 

10 

3.120918 

0.712413 

11 

-0.187885 

-0.225917 

12 

0.090602 

0.764618 

13 

-0.032419 

0.508859 

14 

-0.005370 

0.732557 

15 

0.036480 

0.031661 

16 

-0.890202 

-0.018254 

17 

-0.468842 

0.016737 
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4.  CONTROL  PROGRAM 


A  Fast  Fourier  Transform  (FFT)  subroutine  is  employed 
to  calculate  the  set  tap  coefficients  for  the  equalizer.  The 
equalizer  taps  are  calculated  in  a  recursive  manner. 

The  input  to  the  AMTE  consists  of  a  set  of  N  complex 
correlations  as  a  function  of  frequency.  These  samples  can  be 
expressed  by  the  transfer  function 

H(ju1)  =  K^)  +  j  Q  ( u)  i )  ,  i  =  1  ,2, - ,N 

where  is  the  i  radian  frequency,  I  (w..  )  is  the  in-phase 
frequency  response,  and  Q(ca^)  is  the  quad-phase  frequency  response, 
and  N  is  the  total  number  of  sample  points. 

The  output  of  the  AMTE  is  a  set  of  equalizer  taps. 

For  the  system  at  hand  Nt  is  set  equal  to  17.  These  complex  tap 
coefficients  are  obtained  by  a  truncated  Fourier  series  coeffic¬ 
ients.  The  output  tap  coefficient  will  be  presented  by  its  mag¬ 
nitude  setting  (in  dB's)  and  phase  setting  in  degrees. 

4 . 1  FLOW  CHART 

A  simplified  flow  chart  of  the  main  routine  is  shown 
in  F i gure  4-1 . 

4.2  PROGRAM  IMPLEMENTATION 

The  program  is  implemented  on  the  RADC  HP  2100,  with 
the  operating  system  configured  on  2/6/81.  The  steps  to  be 
taken  for  execution  of  the  program  are: 

After  the  logs  on  the  HP  2100,  the  system  will  return 
the  prompt  sign  At  this  level,  it  is  necessary  to  link 

the  compiled  version  of  the  program  $ Z A : : 4 5 .  To  link,  one  gives 
the  command: 
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FIGURE  4-1.  PROGRAM  FLOW  CKART 
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OF,  AIL  <CR> 


Then  after  computer  responds  with  then  type 

RU,  LOADR,  XX  <CR> 

where  XX  is  the  peripheral  number  of  the  user's 
terminal  (this  may  be  optional)  and  <CR>  is 
the  carriage  return  key 

The  computer  responds  by  giving  a  prompt  of  “/LOADR:''.  After  this 
prompt,  the  user  types: 

RE,  .  ZA  <CR> 

The  computer  then  lists  the  main  line  and  the  subroutines  and 
functions  called  for  by  the  user's  program,  in  this  case,  $ZA::45 
while  l;nking  them  together.  When  the  computer  is  through  linking 
these  together,  it  comes  back  with  another  "/LOADR:".  The  user 
should  then  type 

/E  <CR> 

The  computer  than  links  and  lists  on  the  user's  terminal  as  it  is 
linking,  library  routines  and  functions  necessary  for  executing 
the  user's  program.  The  computer  then  comes  back  with: 

"XX  PAGES  RELOCATED"  "XX  PAGES  REQ'D" 

"NO  PAGES  EMA"  "NO  PAGES  MS EG" 

"/LOADR.-XX  READY  AT  XX  (date,  time)" 

"LOADR-  SEND  and  it  gives  a  prompt 
At  the  prompt,  the  user  should  run  the  program  that  has  been  loaded 
and  linked,  now  called  "AIL:.  The  command  is: 

"RU,  AIL"  <CR> 
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The  computer  then  outputs  anv  messages  and  data  from 
the  program,  as  it  is  running.  The  sequence  of  output  should 
be : 

INITIAL  SETTINGS  AT  STEP  0  SHOULD  BE... 

ENTER  THE  FILE  NAME 

At  this  point,  the  user  types  in  the  name  of  the  input 
file  (containing  a  distorted  sinusoid,  if  created  by  AT10  as  des¬ 
cribed),  for  example  we  could  use  AD000,  one  input  file  from  AT10; 
to  use  AD000  one  types  "AD000"  (return) 

after  the  message  "ENTER  THE  FILE  NAME". 

The  computer  then  prints  out  this  message: 

INTERMEDIATE  SETTINGS  AT  STEP  1  ARE: 

The  computer  then  prints  out  the  attenuator  and  the  line 
stretcher  values  for  the  equalizer,  then  continues. 

If  the  distortion  has  not  been  compensated  in  the  first 
iteration,  the  program  goes  back  to  subroutine  INPUT.  The  program 
repeats  this  loop  10  times,  or  until  the  distortion  is  less  than  a 
threshold  set  to  be  -45  dB,  whichever  comes  first.  When  the  pro¬ 
gram  reaches  the  point  of  exiting  the  loop,  it  prints  out  either: 

"AFTER  10  ITERATIONS  DISTORTION  CANNOT  BE  COMPENSATED" 


or : 


then : 


"DISTORTION  COMPENSATED  AFTER  XX  ITERATIONS" 

"AIL:  STOP" 


The  computer  then  returns  to  the  system  level,  and  the  user  can 
perform  other  tasks,  or  sign  off. 


1 


A  sample  run  of  the  program  implemented  with  input 
data  is  contained  in  Appendix  1. 

4 . 3  COMPUTER  SUBROUTINES 

The  following  discussion  includes  the  entry  points, 
input  arguments  and  output  arguments  of  each  subroutine  in  the 
AMTE  program,  and  what  the  subroutines  or  functions  accomplish. 


1:  Subroutine  INTLZ  ( DMY 1  ,  DMY2 ,  N,  NTAP ,  BETA,  110,  10,  KSW) 

initializes  constants  and  arrays. 

Input:  None 


2  ; 


n 

utput  : 

DMY  1  ,  DMY2 ,  N,  BETA,  10,  KSW 

DMY  1 

A 

Initial  attenuator  (dB)  setting 

DMY  2 

3 

Initial  phase  (degrees)  setting 

N 

Number  of  points  to  be  transformed 

NTAP 

Number  of  taps 

BETA 

Constant  used  in  KASER  subroutine 

10 

Output  device  number 

KSW 

Format  Flag 

KSW  =  0:  the  data  is  expressed  in  decimal  form 
KSW  -  1:  the  data  is  expressed  in  octal  form 

110  Input  device  number  (mag  tape  drive) 

Subroutine  DIAL  (NTAP,  DMY1  ,  DMY2 ,  X,  Y) 

This  subroutine  calculates  attenuator  and  line  stretcher  settings. 
Input:  NTAP,  DMY 1  ,  DMY2 

DMY 1 ,  DMY2  -  updated  attenuation  (dB)  and  phase  (degrees) 
Output:  X,  Y 

X  updated  line  stretcher  (division)  settings 

Y  updated  attenuation  (dB)  settings 
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3:  Subroutine  Head  (X,Y,NTAP,  ISW,  ICNT,  10) 

This  subroutine  displays  attenuator  and  line  stretcher  settings. 
Input:  X,  Y,  NTAP ,  ISW,  ICNT,  10 

X,  Y  line  stretcher  and  attenuation  settings 

NTAP  number  of  taps 

ISW  ISW  =  0;  sidelobe  distortion  <  -45  dB 

ISW  =  1;  at  least  one  sidelobe  >  -45  dB 
ICNT  iteration  counter 

4.  Subroutine  INPUT  (N,  IBUF,  110,  10,  KSW) 

INPUT  defines  input  array  of  in-phase  and  quadrature-phase  data. 

Input:  N,  110,  10,  KSW,  file  name 

File  Name  file  on  which  the  data  is  stored 

Output:  IBUF 

IBUF  an  array  containing  both  I  &  0  data  in  octal  form 

5.  Subroutine  CHNG  (N,  IBUF) 

This  routine  swaps  the  upper  half  of  the  FFT  input  array  with 
the  lower  half. 

Input:  N,  IBUF 

Output:  IBUF 

6.  Subroutine  FORM  (N,  IBUF,  A) 

FORM  will  compute  complex  array  A  =  I  +  jQ  and  also  gain  = 

5  1 o g i q  (I^+Q^)  and  phase  =  tan'^  (Q/I) 

Input:  N,  IBUF 

Output:  A 

7.  Subroutine  INVT  (N,A) 

This  subroutine  computes  the  sequence  to  be  equalized: 
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Input: 

Output 


A.  =  I—  where  i  =  1,2,3 . N 


N  ,  A 


10. 


1 1 


Subroutine  FFT2C  (A.M.IMK)  computes  the  Fourier  transform  of 

the  input  array  A 

Input:  M , A 

Output:  A 

Subroutine  TAPCO  (N,  A,  NTAP ,  COEF,  10) 

TAPCO  truncates  array  A(N)'  into  finite  terms.  The  number  of 
terms  is  defined  as  NTAP  (number  of  taps)  to  specify  an  equal¬ 
izer.  The  variable  10  is  a  device  number  for  the  output  which 
is  normally  computer-dependent.  The  resulting  equalizer  taps 
are  stored  in  COEF  (NTAP). 

Subroutine  KASER  (NTAP,  BETA,  COEF,  10) 

This  subroutine  applies  the  KAISER  window  function  to  the  trunc¬ 
ated  Fourier  transform  output  COEF  (NTAP).  A  functional  sub¬ 
routine  BESSL  (X)  is  required  to  compute  the  window  parameters. 
BESSL  is  the  modified  Bessel  function  of  the  first  kind  with 
zero  order.  The  argument  for  the  BESSL  is  X. 

Subroutine  OUTPT  (0MY1,  DMY2 ,  COEF,  NTAP,  ISW) 

Subroutine  OUTPT  computes  updated  attenuation  (dB)  and  phase 
shift  (degrees) . 

A  DMYLi 

X i  =  20  log1Q  (10*)  ;  A  =  -^r1  +  1 


Y.  =  20  log1Q  ( 1 0D ) ;  B  = 


20 

ATT, 

20 


+  1 


updated  attenuation.  =  20  1  og,  n  (10°) 
X.+Y.  1  lu 

where  C  =  yg  1  +  1 

and  i  =  I  , 2  .....  N 
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DM Y 1  -  previous  attenuation  (dB) 

ATT  -  present  tap  indication  (dB) 

phase  shift  =  phase 

i  =  1  1 

Input:  DMY1  ,  DMY2 ,  C0£F,  NTAP 

DMY1  -  updated  attenuation  setting  (dB) 

DMY2  -  updated  phase  shift  setting  (degrees) 

ISW  -  ISW  =  0;  sidelobe  distortion  <  -45  dB 

ISW  =  1;  at  least  one  sidelobe  >  -45  dB 

12:  Subroutine  DIAL  (NTAP,  DMY1  ,  DMY2 ,  ISW,  ICNT,  10) 

DIAL  was  described  previously. 

13:  Subroutine  HEAD  (X,  Y,  NTAP,  ISW,  ICNT,  10) 

HEAD  was  described  previously. 

14:  Subroutine  END  (ISW,  ICNT,  10)' 

This  subroutine  will  print  a  message  depending  on  what  the 
inputs  are: 

Input:  ISW,  ICNT,  10 

Output:  If  ISW  =  0,  message  is  "Distortion  compensated  after 

so  many  iterations" 

If  ICNT  =  10,  message  is  "After  10  iterations  distor¬ 
tion  cannot  be  compensated" 

15:  Function  ZL0G2:  This  function  calculates  log  to  the  base  10 

of  the  input  quantity  X.  It  takes  the  natural  log  (ALOG)  of  X 
and  multiplies  by  Y,  where  Y  =  1/log  10. 

16:  Function  ZN1 :  This  calculates  arctangent  of  input  quantity 

(in  radians). 
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COMPUTER  PROGRAM  LISTING 


The  following  is  a  listing  of  the  computer  program: 


0  o  0 1 

O  o  he 
0003 
00  04 
0005 
0  0  Ob 
oo  or 
0  0  OS 
0009 
0  0  1  0 
o  o  i i 

0  0 1 3 
0  0 1  3 

0  0 1  4 

0  0 1  5 
0  0 1 6 
0  0 1  7 
00  IS 
0  0  1  9 
0  OS  0 
0021 
0022 
0  02  3 
0024 
0025 
0  02b 
0027 
0028 
0029 
0  030 
0  031 
0032 
0033 
0034 
0035 
00  3  b 
0  037 
0038 
0039 
004  0 
0041 
0042 
004  3 
0044 
0045 


FTN4 

C  REVISED  FILE  l-22-'Sl 

PROGRAM  AIL 

C DUPLEX  A • 25b> > B ( 25b) . CDEF * 25b) 

D 1 MENS I ON  I BUF • 256)  .  DMY 1 ( 1 7 ) . DMY2 <17)«X<17>» Y < 1 ?> 

DIMENSION  IWf (9) 

ICNT=0 

M=8 

C  INITIALIZE  CONSTANTS 

CALL  INTLZ  t DMY 1 »  DMY2  .  N .  NT AP • BETA  *  1 1 □ • IO.KSW) 

C  CALCULATE  INTITIAL  ATTENUATOR  AND  LINE  STRETCHER  SETTING 

CALL  DIAL  (NTAP.  DMY1  *  DMY2.  X.  Y) 

C  DISPLAY  ATTENUATOR  AND  LINE  STRETCHER  SETTINGS 

CALL  HEAD  (X.Y> NTAP» ISW» ICNT. IO> 

DO  10  I  =  1-1 0 
ICNT=  I 

I p ( ICNT. EQ. 2>  GO  TO  20 

C  ENTER  SAMPLE  TRANSFER  FUNCTION  AS  I  AND  &  DATA 

CALL  INPUT  <N.  I  BUF .  1 10.  10.KSW) 

C  SHIFT  DATA  N  2  SAMPLES 

CALL  CHNG  <N» I BUF > 

C  COMPUTE  TRANSFER  FUNCTION  A  (N>  =1  <N>  +JG!  (N) 

CALL  FORM  >  N. I BUF *  A> 

C  COMPUTE  EQUALIZER  TRANSFER'  FUNCTION  l^AcN) 

CALL  INVT i N. A) 

C  PERFORM  FFT  ON  A 

CALL  FFT2C  <A. M. IWk) 

C  TRUNCATE  THE  IMPULSE  RESPONSE  TO  17  TERMS 

CALL  TAPCO  <N. A. NTAP. CDEF  * 10) 

C  PERFORM  WINDOW  FUNCTION  ON  TRUNCATED  SERIES 

CALL  KASEP  (NTAP. BETA. COEF) 

C  CALCULATE  ATTENUATION  (DB)  AND  PHASE  SHIFT  (DEGREES) 

CALL  OUTPT  (DMY 1 . DMY2 . CDEF. NT  AP . I SW) 

C  CALCULATE  ATTENUATION  (TURNS)  AND  LINE  STRETCHER  SETTINGS 

CALL  DIAL  <NTAP. DMY1 . DMY2.X. Y) 

C  DISPLAY  ATTENUATION  AND  LINE  STRETCHER  SETTINGS 

CALL  HEAD  (X. V. NTAP. ISW. ICNT. 10) 

IF 'I SW. EQ. 0) GO  TO  20 
10  CONTINUE 

C  DISPLAY  AFTER  HOW  MANY  ITERATIONS  DISTORTION  WAS 

C  COMPENSATED  .OR  THAT  IT  WAS  NOT  COMPENSATED  AFTER  A 

C  GIVEN  NUMBER  OF  ITERATIONS. 

20  CALL  END -■  ISW.  ICNT.  10) 

9993  fiiR  I TE  (37*  9998 0  > 

99980  F0RMAT(''  AIL  :  STOP") 
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■  *  m 


0047 
0048 
0049 
0050 
0051 
0058 
0053 
0054 
0055 
0056 
0057 
0058 
0059 
0060 
0061 
0068 
0063 
0064 
0065 
0066 
0067 
0  068 
0069 
0070 
0071 
0078 
0073 
0074 
0075 
0076 
0077 
0078 
0079 
0  08  0 
008 1 
0088 
0083 
0084 
0085 


C  CALCUATE  LGG  TO  THE  ERSE  10 

FUNCTION  ZLOG1  <X> 

Y= 0.43489448 19 
IFCX.LE.  0.  >  X=1 .  E-l  0 
ZLOG 1 =Y>ALOG  <X> 

RETURN 

END 

C: 

C  CALCULATE  ARCTANGENT 

C 

FUNCTION  ZN8  <Y.X> 

I F  (  X . EG . 0 . >  X= 1 . E- 1 0 
PI =3. 141598654 
Z=WX 

ZN8=ATAN  vZ'j 
I F  <X . LT . 0>  ZN8-ZN8+P I 
C  WRITE  <37. 976 >  Z.ZN8 

976  FORMAT  <. IX, E18. 5. 1X» El£. 5) 

RETURN 

END 

C 

c 

C  CALCULATE  EESSEL  FUNCTION  OF  THE  FIRST  KIND  WITH  0  ORDER 

FUNCTION  BESSLOO 

Y=X--8 

BELTR=lE-8 
E=  1 . 

DE=  1 . 

DO  1  1=1.85 
DE=DE»Y/FLOAT<:i) 

SDE=DE*DE 

E=E+S0E 

IF  CE^DELTA.  GT.  SDE-1  GOTO  10 
1  CONTINUE 
10  EESSL*E 
RETURN 
END 
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SWAP  THE  UPPER  HALF  OF  THE  ARRAY  WITH  THE  LOWER  HALF 
SUBROUTINE  CHNG  CN* IBUF> 

INTEGER  IBUF  Cl  >  * IE'UFl  (30 CO 
K.  1  =  CN+D 
K2=K1-1 
k>N-'£ 

BO  10  I =K 1 » N 
J-I-K 

IEUF1 <J>=IBUFCI> 

CONTINUE 
DO  £0  1*1 ,K 
J=I+K1 

IBUF1 <J>  =  IBUF  <I> 

CONTINUE 
BO  30  1  =  1  *N 
IBUF CI>  =IBUF1 <I> 

CONTINUE 

FORMAT  C1X, 1£CI6» IX) > 

ZX=1 

CALL  TYPE  CZX) 

RETURN 

END 


I 


0110 

C 

01 1 1 

c 

0118 

c 

0113 

c 

0114 

c 

01 15 

0116 

0117 

0 1 1 8 

01 19 

0180 

0 1 8 1 

0188 

0133 

0134 

0135 

0186 

0137 

0 1 88 

0139 

0130 

0131 

0133 

0133 

0134 

0 1 35 

0136 

0137 

0138 

0139 

0140 

0141 

0148 

1  0 

0143 

c 

0144 

c 

0145 

979 

0146 

0147 

0148 

0149 

PERFORM  POLYNOMIAL  FUNCTION  TD  CONVERT  FROM  ATTENUATION  <DB  ' 
TO  ATTENUATION  'DIAL  TURNS)  AND  FROM  PHASE  SHIFT 'DECREES > 
TO  LINE  STRETCHER  SETTINGS  'DIVISIONS) 

SUBROUTINE  DIAL  CNTAP  >  DMY1 * DMY2  *  X*  Y> 

D I  MENS- 1  ON  A  0  ( 1 7  :■  *  A 1  C 1 7)  *  Ac!  <  1 7)  *  A3  (.17)  * 

1  DMY Kl).  DMY8  C  1  >  *  X  cl)  *  V  C 1  > 

DAT  A  AO'  -  •  3 1  03  *  — .  7  c  6 1  *  — .  935 1  *  — .  688 1  *  — .  86  07 * —  1 .  £897 * 

1  —  1 . 1816*—  1 . 1316*  0 . *  - . 8644 *  — . 84 16*  —  1 .  038 1 » — . 91 69 » 

3-8 .  1413»~.  84c 8  *  -  1  •  034c' *  — .  7 3 4 3 
DATA  A 1 /— . 1 166?  —  .  07388*-.  1116*-. 06838 *  - . 1 086 *  - . 1 8938 * 

1-.  1  075* -.  1  075*  0.  *-.  08305* 07806* -.  1176*-.  0946*-.  1399* 

£-. 05173*-. 1193*-. 0777 1 / 

DATA  A8v-. 0009174* . 001395*-. 0003383* . 001839* . 0003681  * 

1 . 0  0  0  071 63*-.  0  0 03634 *  - .  0 0 0  3684 *  0 . * .  0  0 03  096  *  - .  0  0 08863  * 

3-.  00003543*  .  000673* .  0006748*  .  008387*  -.  001 177*  .  001385'- 
DATA  A3 • '  - .  0  0  0  0 3 8 3797* .  0  0  0  0  0377*-.  0  0  0  0 8436*.  0  0  0  0 1 1 0 7 5  * 

1  - .  0  0  0  0318*-.  0  0  0  0  3 3  3 18*- .,  0  0  0  08773*-.  0  0  0  0 8773*0.*-.  0  0  0  0 8  5  0  9 6  * 
£-. 00003957*-. 0000373*-. 000016696*-. 00001571* . 00001876* 

3- .  0  0  0  04344  *  - .  0  0  0  0  083 1  €•/ 

DO  10  1=1 *NTAP 
J=I 

IF  <1 . LE. 8) J=9-I 

Y 1  =  A  0  C  I)  +A 1  ( I  >  ♦DMY  1  <  J>  +A8  Cl)  ♦DMY1  (.J)*DMY1  CJ) 

Y3-A3  Cl) ♦DMY 1  CJ) ♦DMY1 CJ) ♦DMY 1 CJ) 

Y  C I ) =Y1 +Y8 

X  Cl)  =DMY£  CJ)  ♦3.634''  <360.  ♦!.  091) 

I F  CY  C I )  .  LT’.  0 .  0)  Y  C  I  )  =  0 .  0 
I F  CY  C 1 )  .  6T .  1  0 .  0)  Y  C I )  =  1  0 .  0 
IF  CX  C I )  .  LT .  0 .  0)  X  C I )  =  0 .  0 

IF  CXCD.6T.  11.  0)  X  CD  =1 1.0 
CONTINUE 

MR I TE  C37 *  979)  CX Cl) . 1  =  1* NTAP) 

WRITE  ( 37*979)  (Y (I) *  1  =  1*  NTAP) 

FORMAT  < 1 X*  4  CE 1 3. 5* 1X*E13.5* IX) > 

2X=£ 

CALL  TYPE  CZX) 

RETURN 

END 
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0 1 5  0  c 

|.'|  151  C 

0152  C  PRINT  DUT  WHETHER  DISTORTION  IS  COMPENSATED  AFTER  THE 

0153  C  LAST  ITERATION 

0 1 54  SUBROUT  I NE  END  < I S W. I C NT . 1 0> 

0 155  IF  <1 SW .  EG' .  0>  GO  TO  10 

0 1 56  WR I TE  <  3  7 .  £  0;-  I NC T 

0157  £0  FORMAT  < ✓ , 3X  »  5HAFTER » 1 3 »  3 1 H  ITERATIONS  DISTORTION  CAN  NOT. 

0158  114HEE  COMPENSATED.--'. 

0159  GO  TO  30 

0160  10  CONTINUE 

0161  ICNT=ICNT-1 

0 1 6£  WR I TE (37. 4  0> I C  NT 

0163  40  FORMAT!.  .3X.28HD I STORT I ON  COMPENSATED  AFTER. I£. 1 1HITERAT ION 

0164  30  CONTINUE 

0165  ZX=3 

0166  CALL  TYPE  <!ZX> 

0167  RETURN 

0168  END 

0169  C 
0170  C 
0171  C 

0173  SUBROUTINE  FORM <N. IBUF. A) 

0173  COMPLEX  A<1> 

0174  DIMENSION  IBUF<1> 

0175  C  COMPUTE  COMPLEX  ARRAY  A=I+JQ 

0176  DO  10  1  =  1. N 

0177  I A1 =  1 EUF  < I > /£56 

0178  I B 1  =  I BUF  < I > -£56* I A 1 

0179  IFdAl.GE.  1£8>  IAl  =  IAl-£56 

0180  IFdBl.GE.  1 28>  IBl  =  IBl-£56 

0181  XI “FLOAT  < IA1> ' 128. 

0182  XQ=FL0AT<IB1>.'1£8. 

0183  A  < 1 > =CMPLX  <X I . XG) 

0184  10  CONTINUE 

0185  998  FORMAT  < IX. 4 (El£. 5. 1X.E12. 5. IX) > 

0186  ZX=4 

0187  CALL  TYPE  <2X> 

0188  RETURN 

0189  END 

(.  1  «•  0  C 
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T 


0191 

0 1 9  c' 

c 

SUBROUTINE  HERD  <X*  Y»  NTAP» I SW. I C  NT « 1 0> 

0 1 93 

n i mens i on  x a>»  v  <: i  • « zr  <: l ?> 

0194 

DRTR  ZA-4H1  H= .  4H3  H=.4H3  A=.4H4  A=«4H5  A=j4H6  A=i>4H7  R=» 

01 95 

1 4H8  R=  *  4H  M  =  *4H  1  B=*4H3  B=-4H3  E=. 4H4  E-.4H5  E-.4H6  B=* 

0 1 96 

1 4H7  B= . 4H8  B=" 

0197 

WRITE '37. 145> 

0 1 96 

145 

FORMAT * 1H1 > 

0 1 99 

I  CNT  1  =  1  CNT  +•  1 

03  00 

bO  TO  *.  1  0 .  3 0 .  3 0 *  3 0 » 3 0 » 3 0 *  3 0 *  3 0 »  3 0  j  3 0 *  30'*  ICNT1 

03  01 

1  0 

CONTINUE 

03  03 

WRITE • 37- 1 50> I CNT 

03  03 

j” 

DISPLAY  ATTENUATION  AND  LINE  STRETCHER  SETTINGS 

0304 

1  5  0 

FORMAT*  18X.35HINITIAL  SETTINGS  AT  STEP  »I1»11H  SHOULD  EE: 

03  05 

GO  TO  305 

03  06 

30 

CONTINUE 

0307 

I SW 1 = I SW+ 1 

03  06 

bO  TO  *  ij .  c‘  U  0  ■’  I  S W 1 

03  09 

3  0  0 

CONTINUE 

03 1  0 

WRITE <37. 155> I CNT 

03 1 1 

155 

FORMAT  •:>  ,  18X.39H INTERMEDIATE  SETTINGS  AT  STEP»I£»5h  ARE:  > 

0313 

GO  TO  305 

03 1  3 

3  0 

CONTINUE 

03 1  4 

WRITE  *37. 1 60> I CNT 

03 1 5 

1 6  0 

FORMAT*/. SIX. 33HFINAL  SETTINGS  AT  STEP  » 12, 5H  ARE:  > 

03 1 6 

305 

CONTINUE 

03  1  7 

WRITE *37. 965  > 

03 1 8 

965 

FORMAT  *  -v  1  OX.  1  7HATTENUATOR  (TURNS )  »31X»  19HLINE  STRETCHER  (.DIV 

03 1  9 

DO  310  1  =  1.  NTAF' 

0330 

WR I TE  *  37 .  17 ZA  <  I )  » Y  <  I )  •  ZA  < . I •  X  *.  I } 

0331 

170 

FORMAT  < 1 0X«  H4»  F6. 3»  34X . A4»  F4. 1> 

0333 

310 

CONTINUE 

0333 

ZX=5 

0334 

CALL  TYPE  CZX) 

0335 

RETURN 

0336 

0337 

C 

END 
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0228  C 

0228  SUBROUTINE  INPUT  <N» IBUF. I IO» IO*kSW> 

0230  D I MENS I  ON  I BCB <256> « NAM <3> # I BUF c £58) 

023 1  D I  MENS I ON  I BUF 0 < 1 £8> *  I BUF 1 < 1 £8> 

0232  WRITE  <37* 1 0> 

0233  10  FORMAT  <>'•  '  ENTER  THE  FILE  NAME 

0234  REND  < 37 . 15> NAM 

0235  15  FORMAT  <.3A£> 

0238  I L=M 

0237  CALL  OPEN  <  I  DC B.  IERR.  NAM> 

0238  IF (I ERR. LT. 0>  GO  TO  900 

0239  C  ENTER  DATA  IN  OCTAL  FORMAT 

0240  CALL  READFCIDCB. IERR. IBUF0> 

024 1  CALL  RSfiliF  < I  BCB  .  I  ERR  .  I  BUF  1  > 

0242  IF  (IERR. LT. 0>  GO  TO  910. 

0243  CALL  CLOSE CIDCE. IERR) 

0244  C  FILL  ARRAY  WITH  IN-PHASE  AND  QUADRATURE-PHASE  COMPONENT 

0245  DO  80  J=  1 »  1 28 

0248  I BUF  < J> = I BUF  0  <  J> 

0247  80  CONTINUE 

0248  DO  70  K=1 .  128 

0249  IEUF <K>1£8)=IBUF1 <K> 

0250  70  CONTINUE 

0251  998  FORMAT  <  1 X  >  1  £  <  1 8 . 1 X)  > 

0252  60  TO  999 

0253  900  WRITE <37. 30> IERR 

0254  30  FORMAT  <>» "  FMP  ERROR  "14) 

0255  GO  TO  999 

0258  910  WRITE <37. £0> 

0257  £0  FORMAT  '  ERROR  READING  THE  FILE  O 

0258  999  ZX=8 

0259  CALL  TYPE  <ZX> 

0260  RETURN 

028 1  END 

0262  C 
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0263  C 
0264  C 
0265 

02  6  b 
02'  6  7 
0268 
0269 
027  0 
0271 
0272 
0273 
0274 
0275 
0276 
0277 

0278  1 0 

0279 
0280 
0281 
0282  C 
0283  C 
0284 
0285 
0286 
0287 
0288 
0289  C 


INITIALIZE  CONSTANTS  AND  ARRAY 
SUBROUT I NE  I NTLZ t DMY 1 *  DMY2 * N * NTAP «  BETA* 1 10? KSW) 

D I  MENS I ON  DMY 1 < 1 >  * DMY2 Cl >  ,  A  Cl 7> *  B  C 1 7> 

DATA  A.  -59. 964*  -59. 964  *  -51 . 072* -56. 743  *  -77.  661  >  -6.0. 55 1  * 
1  —  72 ■  7I4<  57,  Z‘ 5  o'  *  0 •  '  —  5 8 .449*  — 6c' .  143*  —  55.  o'  s  8  *  5 »  ■  r  4 1  * 

2-54 .419*  -69 . 5  04 *  —57 .  047  *  —58 . 965/ 

DATA  B/659 .  23 *  756 . 49 >  8 1  0 . 53  *  983 .  44  *  162.  1*399.86* 

167  0 .  03 » 859 .  16*443.  09  *  1  c'9 .  66  *  399 .  86.  •  56- 1 . 96.  *  778 .  1  .* 

254 .  04  *  151.3*345. 82  *  670 .  03/ 

N=256 
NThP=1  7 

DO  10  1=1* NTAP 
DM VI C I >  =A  C I > 

DMY2 CI> =  BCI> 

CONTINUE 
BETA=1. 0 
1 10=37 
10=37 

FOR  DECIMAL  I«<Q*KSW»0 
FOR  OCTAL  l&Q* KSW-1 
KSW=1 
ZX=7 

CALL  TYPE  <ZX> 

RETURN 

END 
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1 


0890 

C 

089 1 

C 

CALCULATE  THE  RECIPROCAL  OF  THE  TRANSFER  FUNCTION 

0292 

SUBROUTINE  INVT(N.R) 

0293 

COMPLEX  A<1> 

0294 

DO  10  I  =  1  •  N 

0295 

A  f  I  >  » 1 .  /  (flU)  ♦FLOAT  (N>  > 

0296 

1  0 

CONTINUE 

029? 

c 

MR  I  TE  (37 «  998)  (A  ( I )  .  I  =  1  > N.) 

0292 

FORMAT  1 1 X .  4  <  E 1 8 . 5 . 1X.E18.5. 1X>> 

0299 

ZX=8 

0300 

CALL  TYPE  <ZX) 

0301 

RETURN 

03  02 

END 

0303 

c 

0304 

c 

0305 

c 

APPLY  KAISER  WINDOW  FUNCTION  TO  SMOOTH  OUT  THE  RIPPLE 

030c 

c 

TO  TRUNCATION  OF  THE  FOURIER  SERIES 

0307 

SUBROUTINE  KASER (N.B.C) 

0308 

COMPLEX  C  < 1 > 

0309 

IF  CB.EQ.  0.0)60  TO  9999 

03 1  0 

NVS*  (N+l )  ‘'2 

031 1 

DO  10  1  =  1  >N 

0318 

Z8=FL0AT (I-NVS) AFLOAT  <N'8> 

0313 

Z3=B*SQRT  <1 . -Z2*Z2 ) 

0314 

C  (I >  =C  ( I >  ♦BESSL  <Z3>  BESSL  (B> 

0315 

1  0 

CONTINUE 

0316 

c 

WRITE  (37.979)  (C(D.I  =  1.N> 

03 1 7 

979 

FORMAT  (  1 X . 4 ( 1 X . E 1 8 . 5 . 1 X . E 1 8 . 5> ) 

0318 

9999 

ZX=9 

0319 

CALL  TYPE (ZX) 

0380 

RETURN 

0381 

END 

0388 

C 
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0323 

0324 

0325 

0326 

0327 

0322 

0322 

033  0 

0331 

0332 

0333 

0334 

0335 

03  36 

0337 

0338 

0  339 

034  0 

0341 

034  2 

034  2 

0344 

0345 

0346 

0347 

0348 

0349 

0350 

035 1 

0352 

0353 

0354 

0355 

0356 

0357 

0358 

0359 

0360 

0361 

0362 

0363 

0364 

0365 

0366 

0367 

0368 

0369 

0370 

0371 

0372 

0373 


C 

.  UBROUT  I  HE  OUTPT  ' DMY 1 ♦ DMY2« C OEF  •  NTRP.  1 1  Id;. 

COMPLEX  C  GEF  •:  1  > 

D I MENS  ION  DMY 1 <1>. DMY2 <l>»ftTT<17), PH 3 il7) 

DO  10  1=1. NTRP 
XR=REhL  *COEF  •:  I  )  > 

XI  =R  I  MR6  <  C  OEF  >  I  :•  > 

C  COMPUTE  RTTENURTION  OF  DISTORTION 

RR=XR^XR+XI ♦XI 
RTT  <I>  =1 0. ♦ZLC61  <RR> 

C  COMPUTE  THE  PHRSE  SHIFT  OF  TRP  COEFFICIENTS 

PHS  <  I  >  =.:N2  <  XI  .  XP>  ♦45.  -'PtTRN  Cl .  > 

ZZ=ZN2<XI • XR • 

C  iitR  I  TE  <27»978>  PHS  <  I  >  *  ZZ 

978  FOF'MRT  •  1 X  •  E  1 2 .  5 .  1 X  ••  E 1 2 . 5> 

10  CONTINUE 

DO  20  I = 1 • NTRP 

RTT  i  I  =RTT  >  I  >  -RTT  <9> 

PHS i I > =PH: < I ' -PHS < 9> 

80  CONTINUE 
I  S  1.1=0 

NT  RP 1 =NT  RP- 1 
DO  30  I=1.NTRP1 
J=I 

IFCI.6E.9)  J=I  +  1 
IF  (RTT  <J>  .  LE.  -45.  >  GO  TO  15 
I  SW=  l 
60  TO  30 
15  CONTINUE 
I SU= I SW 
30  CONTINUE 

I F  1 1 S W .EC1.  0 >60  TO  4 0 
DO  50  1=1 -NTRP 

X=20. ♦ZL06 1  <  < 1 0. ♦♦ ' DMY 1 < I>  2  0. >  >  +1> 

Y=2 0  ♦ZL06 1 < < 1 0 . ♦♦ <  RTT c I > /£ 0 . > >  + 1 > 

DMYi  <  I  >  -20.  ♦ZLOG 1  '■  k  1  0.  ♦♦  *  <X+Y>  •••'20.  >  >  — 1> 

DMY 2 C I > =DMY2 < I > +PHS < I > 

50  CONTINUE 
40  CONTINUE 

WRITE  c37» 996 > 

996  FORMRT  <1X«  1 4HRTTENURT0R <D£0  >37X»  19HLINE  STRETCHER  <  DE6>  ■>  -■  > 
*  888  ZX= 1 0 

DO  881  1  =  1  * NTRP 

UR  I TE  <37 • 977>  DMY 1 < I > • DMY 2  ( I > 

977  FORMRT  < 1 X. E12. 5< 37X. El  2. 5> 

881  CONTINUE 

979  FORMRT  < IX. 4 <J X» E12. 5* IX. E12. 5> > 

CRLL  TYPE  <ZX> 

RETURN 

END 

C 
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0374  C 

0375  SUBROUTINE  TAPCO  (N» A»  NTAP*  COEF »  IO> 

0376  COMPLEX  A  v  1  >  « COEF < 1 > 

0377  NV 3=  CNTRP+ 1 > <  3 

0378  NV31 =N+NV3+1 

0379  NV1=NV3+1 

0380  C  TRUNCATE  THE  ARRAY  A  <N>  INTO  A  GIVEN  NUMBER CNTAP)  OF 

0381  C  TERMS 

0383  DO  30  1  =  1  * NV3 

0383  J=NV£-I+1 

0384  COEF<I>=AOJ> 

0385  30  CONTINUE 

0386  C  OBTAIN  THE  UPPER  HALF  OF  THE  COEFFICIENTS  CENTERED 

0387  C  AROUND  <NTRP+1> ✓£ 

0388  DO  30  I =NV 1 «  NTRP 

0389  J=NV3 1 - I 

0390  COEF<I>-A<J> 

0391  30  CONTINUE 

0393  C  WRITE  <37»979>  CCOEF <I> » 1=1 » MTAP) 

0393  979  FORMAT  C1X> 4 <1X» E13. 5» 1X> E13. 5>  ) 

0394  ZX=1 1 

0395  CALL  TYPE  <ZX> 

0396  RETURN 

0397  END 
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0398 

0399 

0400 

0401 

04  08 

0403 

04  04 

0405 

04  06 

040? 

0408 

0409 

04 1  0 

0411 

04 1  £ 

0413 

0414 

04 1 5 

04  1 6 

04 1 7 

04 1 8 

0419 

04  c  0 

04c  1 

0438 

04  £3 

04  £4 

04£5 

04£6 

04  £? 

04£8 

04£9 

0430 

0431 

043£ 

0433 

0434 

0435 

0436 

043? 

0438 

0439 

0440 

0441 

04 4 £ 

0443 

0444 


C 

C 

SUBROUTINE  FFT£C  *A*M*IWK> 

INTEGER  M*IWKCl) 

COMPLEX  H  <1  ;■ 

C 

C 

INTEGER  I *  ISP»  J»  JJ»  JSP •  K  » K0»  K 1 » K£ » K3 *  KB » KN » MK » MM*  MP*  N» 
1  N4 • N8 *  N£  * LM «  NN  *  Jk 

R'EAL  RAB*  C  1 *  Cc'j  Co?  SI  •  S£«  S3*  CK *  3K  >  SO?  A0  *  A1 *  A£>  A3* 

1  B0* B 1  * B£* B  3 * TWO PI • TEMP* 

C  £  ZERO*  ONE*  Z0  c£>  *  Z1  C£>  *Z3  C£> 

£  ZERO  *  ONE  •  Z  0  C  £ )  *  Z 1  <  £>  *  Z3  C£>  •  Z£  *  £> 

COMPLEX  ZAO*  ZA1  >  ZA£«  ZA3*  Ak£ 

EQU I  VALENCE  CZA0  *  Z0  <  1  >  )  *  CZR  1  *  Z  l  <  1  >  >  *  CZA£  *  Z2  C  1 )  >  * 

1  CZA3* Z3 < 1 ) > > CR0*  ZO  Cl>  >  * <B0»  ZO  <£>  >  » Cfil » 21 <l>  >  » 

£  <  B 1  *  Z 1  <  £)  )  *  <  H£*  Z£  C 1 >  «  C££*  2£  C£)  >  »  c A3  *  Z3  c 1 ) ) * 

3  (B3* Z3 <£> > 

BATA  SO'  .  7071  06781' 

BATA  S'  K  / .  3  8  £ 6 8 3 4  3 £  / 

DATA  C  K  .  9  £ 3 8795 S;  / 

BATA  TWQP I -'6.  £83 1 853 1  / 

BAT  A  ZERO/ 0 • 0/ *  ONE/ 1 . 0/ 

MP  *  M+l 
N  =  £♦♦»<! 

I WK  Cl)  =  1 
MM= CM/£) ♦£ 

K  N  =  N+l 
BO  5  I  =£  •  MP 

I  WK  ci>  = i  wk  <  i  - 1 )  + 1  uik  c  i  - 1  > 

5  CONTINUE 

PAD  =  TWOPI/N 
MK  =  M  -  4 
KB  =  1 

IF  CMM  .EG.  M>  GO  TO  15 
K£  =  KN 

K0  =  IWKCMM+1)  +  KB 
1 0  K£  =  K£  -  1 
K0  =  K0  -  1 
AK£  =  A  cK2> 

A  <K£>  =  A  CKO)  ~  AK£ 

A  CKO)  =  A  CKO)  +  AK£ 

IF  CKO  .GT.  KB)  GO  TO  1 0 
15  C 1 =ONE 

SI  =  ZERO 
JJ  *  0 
K  =  MM  -  1 


4-20 


0445 

J  =  4 

0446 

IF  (K  .GE.  1>  GO  TO  30 

0447 

GO  TO  70 

0448 

£0 

IF  (IWK<J>  -GT.  JJ>  GO  TO  25 

0449 

JJ  =  JJ  -  IWK  (J> 

0450 

J  =  J  -  1 

0451 

IF  (IWK(J>  .GT.  JJ>  GOTO 

£5 

0452 

JJ=JJ-IWK  CJ> 

0453 

J=J-1 

0454 

k  »  K  +  2 

0455 

GO  TO  20 

0456 

25 

JJ  =  IWK (J)  +  JJ 

0457 

J  =  4 

0458 

30 

ISP  =  IWK(K> 

0459 

IF  (JJ  ,EQ.  0>  GO  TO  40 

0460 

C2  =  JJ  ♦  ISP  ♦  RRD 

0461 

Cl  =  COS  <C2> 

0462 

SI  =  SIM (C£> 

046? 

35 

C£  =  Cl  ♦  Cl  -SI  ♦  SI 

0464 

S£  *  Cl  ♦(SI  +S1> 

0465 

63  =  C£  ♦  Cl  -S2  »S1 

0466 

S3-  C£  ♦  SI  ♦  S2  ♦Cl 

0467 

4  0 

JSP  *  ISP  ♦  KB 

0468 

C 

WRITE  (37  r 997>  S1»S£?S3 

0469 

C 

WRITE  (37 > 997>  C1»C£»C3 

047  0 

DO  50  I  =  1  »  ISP 

0471 

KO  =  JSP  -  I 

0472 

K1  =  KO  ♦  ISP 

0473 

K£  =  K1  +  ISP 

0474 

K3  -  K2  +  ISP 

0475 

ZRO  *  R  (K0> 

0476 

ZR1  =  fi(Kl> 

0477 

ZR£  *  R (K£> 

0478 

ZR3  =  R (K3) 

0479 

IF  (SI  .EQ.  ZERO>  GO  TO 

45 

0480 

TEMP  *  R1 
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0481 

81  =  ft 1  ♦  Cl  -B1  ♦SI 

0482 

El  =  TEMP  ♦  SI  +  B1  ♦  Cl 

0483 

TEMP  =  82 

0484 

82  =  82  ♦  C2  -  B2  ♦  S2 

0485 

B2  =  TEMP  ♦  S2  +  B2  ♦  C2 

0486 

TEMP  =  83 

0487 

83  =  83  ♦  C3  -  B3  *S3 

0488 

B3  =  TEMP  ♦  S3  +  B3  ♦  C3 

0489 

45 

TEMP  =  80  +  82 

0490 

82  =  80  -  82 

0491 

80  =  TEMP 

0492 

TEMP  =81  +83 

0493 

83  =  81  -83 

0494 

81  =  TEMP 

0495 

TEMP  =  BO  +  BE 

0496 

B£  =  BO  -  B2 

0497 

BO  =  TEMP 

0498 

TEMP  =  El  +  B3 

0499 

B3  =  B1  -  E3 

0500 

B1  =  TEMP 

0501 

8Z 1=80+81 

0502 

8Z2=B0+B1 

0503 

8Z3=80-81 

0504 

8Z4=B0-B1 

0505 

8Z5-82-B3 

0506 

8Z6=B2+83 

0507 

8Z7=8£+B3 

0508 

8Z8=B£-83 

0509 

997 

F0RM8T  aX.4<E12.5.  1X>> 

05 1  0 

RCK0>  =  CMPLX  C8Z1 » 8Z£> 

0511 

8<K1>  =  CMPLX C8Z3.8Z4> 

0512 

8  (K2>  =  CMPLX C8Z5. 8Z6> 

0513 

8  <K3>  =  CMPLX <8Z7 . 8Z8) 

0514 

999 

FORM8T  ClX.'R  COEFFICIENTS’') 

0515 

50 

CONTINUE 

0516 

C 

WRITE  <37.  998>  8<K0>  »8<K1>  »8<K2>  »8<K3> 

0517 

C 

WRITE  <:37» 995>  K0»K1»K£»K3 

0518 

995 

F0RM8T  aX.4CI4-.lX>> 

0519 

998 

F0RM8T C 1 X. 4 <E 12. 5. IX. El 2. 5. 1X> > 

0520 

IF  <K  .LE.  1>  60  TO  55 
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*-1  OJ  00  If  U"|  Mj  f'-  CO  O'  *1  OJ  CO  't  IT*  •£>  t-  CC'  O'  O  «-i  CU  0?  ’t  BO  •£)  N  CO  (Ji  C>  «  OJ  00 

cu  oj  oj  oj  oj  oj  oj  oj  oj  oo  co  co  co  co  co  oo  oo  co  oo  ■>*  tj-  rf  >t  -t  tj-  tj-  <t  tj-  ir>  ifi  no  in 

no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no  no 


k  =  K  -  £ 

GO  TO  30 

55  KB  =  K3  +  ISP 

IF  CKN  .LE.  KB>  GD  TO  70 
IF  <J  .  NE.  1>  GO  TO  60 
K  =  3 
J  =  MK 
GO  TO  £0 
60  J  =  J  -  1 
C£  =  Cl 

IF  <J  .NE.  £>  GO  TO  65 
Cl  =  Cl  ♦  CK  ♦  SI  ♦  SK 
SI  =  SI  ♦  CK  -  C£  ♦SK 
GO  TO  35 

65  Cl  =  CC1-S1>  ♦  SG» 

SI  =  CC£  ♦  SI)  ♦  SQ 
GO  TO  35 
70  CONTINUE 

C:  WRITE  <3 7?  998 >  CfiCKZ>  >KZ=1  »£56> 

IF  CM  .LE.  1>  GO  TO  9005 

MP  =  M  +  1 
JJ  =  1 
I WK  Cl)  =  1 
HO  75  I  =  £*  MP 

I WK  CI>  =  I WK  C I - 1 > ♦£ 

75  CONT I NUE 

N4  =  I WK  <MP-£> 

IF  CM  .GT.  £)  N8  =  IWKCMP-3) 

N£  =  I WK  CMP— 1 > 

LM  =  N£ 

NN  =  I WK  CMP)  ♦  1 
MP  =  MP  -  4 
J  -  2 
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.  .  MmiitUf  SJHm*  . 


.  .  Vi 


0354 

80 

JK  =JJ  +  N2 

0555 

ftK2  =  ft  <J> 

0556 

ft  C  J)  =  ftCJK) 

0557 

ft  C  JK)  =  ftK£ 

0558 

J  =  J+l 

0559 

IF  CJJ  .  GT.  N4)  GO  TO  85 

056  0 

JJ  =  JJ  +  M4 

0561 

GO  TO  105 

0568 

85 

JJ  =  JJ  -  N4 

0563 

IF  CJJ  .GT.  NS)  GO  TO  90 

0564 

JJ  =  JJ  +  N8 

0565 

GO  TO  105 

0566 

90 

JJ  =  JJ  -  N8 

0567 

k  =  up 

0568 

95 

IF  CIWKCIO  . 6E.  JJ>  GO  TO 

1  00 

0569 

JJ  =  JJ  -  IWKCK) 

0570 

K.  =  K  -  1 

0571 

GO  TO  95 

0572 

1  0  0 

JJ  =  IWK  CK>  +  JJ 

0573 

105 

IF  CJJ  .LE.  J>  GO  TO  110 

0574 

k  =  Nft  -  J 

0575 

JK  s  Mi  ~  JJ' 

0576 

ftK£  =  ft  C  J> 

0577 

ft  C J)  =  ft  CJJ) 

0578 

ft  CJJ)  *  ftK2 

0579 

ftK2  *  ft  CK> 

0580 

ft  CIO  =  ft  C  JK) 

0581 

ft  C JK)  =  ftK2 

0588 

110 

J  »  J  ♦  1 

0583 

IF  CJ  .LE.  LM)  GO  TO  80 

0584 

9005 

ZX=  1 2 

0585 

CftLL  TYPE  CZX) 

0586 

C 

WRITE  C37 • 998)  Cft  CKZ) >  KZ«1 

..  256 

0587 

RETURN 

0588 

END 

0589 

C 

0590 

C 

0591 

SUBROUT I NE  TYPE  CZX) 

0592 

C 

WRITE  C37 • 20) ZX 

0593 

20 

FORMftT </i '  SUBROUTINE  ENTERED 

0594 

RETURN 

0595 

END 

:  \ 
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4.5 


COMPUTER  SOFTWARE  VERIFICATION 


To  assure  that  the  software  works  properly,  the  follow¬ 
ing  verification  procedure  has  been  established.  Before  running 
the  program  on  real  data,  it  is  recommended  that  the  program  be 
run  on  an  idealized,  sinusoidal  waveform  with  slight  distortion 
added . 

Two  cases  of  distortion  will  be  introduced;  the  first 
one  involves  amplitude  only,  and  the  second  one  involves  phase 
only.  The  distortions  are  presented  as  sinusoids  of  various 
amplitudes  but  always  small  enough  so  as  to  permit  the  use  of 
the  paired  echo  theory  concepts  to  be  valid.  Also,  the  period 
of  the  sinusoid  is  always  harmonically  related  tr>  the  bandwidth 
of  the  equalizer  network  and  as  such  is  consistent  with  the  con¬ 
cept  of  representing  the  distortion  as  a  Fourier  series. 

The  test  cases  of  distortion  are  separated  into  two 
groups,  one  of  which  assumes  only  amplitude  distortion  and  the 
other  of  which  assumes  only  phase  distortion,  but  of  varying 
levels.  The  FFT  presents  the  coefficients  for  the  taps  which 
in  effect  represent  the  amplitude  and  phase  of  the  echo  associ¬ 
ated  with  each  tap.  These  coefficients  are  then  compared  to  the 
calculated  values  based  on  the  paired  echo  theory. 

CASE  1:  A  sinusoidal  amplitude  distortion  of  various  peak  to 
peak  values. 

| G|  =  10  log1Q  v/ I2+Q2  =  5.1og10  (I2+Q2)  =  amplitude 
0  =  tan'1  ^  *  phase 
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_ III 


To  simplify,  we  consider  the  effects  of  only  amplitude  distortions 
so  that  assuming  no  phase  distortion  and  letting  Q  =  0  we  get 

|G|  =  5  log1Q  I2  =  10  log1Q  I 

0  =  tan-1  (|)  =  0 

Since  a  sinusoidal  distortion  is  needed: 

-1  n-  =  10  log1Q  I 

where  AX  =  incremental  frequency 

N  =  number  of  points 

K  =  scale  factor  (for  K  =  1,  numerator  goes 
0  to  -2  or  2  dB 

and  10  log  I  represents  the  distortion  level  in  dB  ( peak-to-peak )  . 

In  order  to  fit  the  RADC  data  form  so  that  -128  <  I  <  128,  we  multiply 
the  given  values  of  I  by  128. 

TEST  1:  Amplitude  distortion,  0.5  dB  peak-to-peak  with  first 
echo  pair 

-0.5  =  L^I  =  io  log1Q  therefore  K  =  4 

For  215  points  in  2n  radians  or  360°,  there  AX  =  and  0  <N  <214* 

i  =  io  [s-la  ffigj ~-l] 

From  the  test  run  we  get  the  first  echo  peak  value  distortion  to  be 
-30.44  dB.  Analytically  (for  small  amp  distortion)  we  have: 

♦Note  that  the  number  of  sample  points  is  not  germane  to  the  analysis 
(as  long  as  it  is  a  large  number)  and  relates  only  to  the  expected 
number  of  samples  to  be  used.  We  later  change  this  to  256  points. 
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®  1 

20  IoSiq  (’+  "— )  =  amP^itucie  nipple  in  dB  (peak  value) 


a,  1  a1 

for*  ~  <<  1  ;  j  -1 -  s  peak  echo  amplitude 
ao  c  ao 


20  1o9iq  )  *  echo  in  dB 
o 


20  log1Q  (1+  j-L)  -  0.5  dB 


ai  lo 

,  jo - | O -  =  0.0296 


echo  *  20  log1Q  (jj-)  =  -30.566  dB 


TEST  2:  Distortion  is  introduced  in  the  second  echo  pair  with  an 
amplitude  ripple  of  0.5  dB 


I  ■  10 


where  Ax  ■ 


2x360 

“2T5- 


and  0  <  N  <  215 


From  the  test  run  we  get  second  echo  peak  value  to  be  -30.53  dB  where 
the  calculated  value  is  the  same  as  before  -30.566  dB. 

TEST  3:  Distortion  is  introduced  in  the  first  and  second  echo  pairs 
with  an  amplitude  ripple  of  0.5  dB 


I  *  (10  [- 


sin  (Ax,N)-l 


-3  +  10  [- 


sin  (Ax-N)-l 


where  Ax]  -  |y|;  A x2  =  2^-|°  and  0  <  N  <  214 
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Because  we  are  summing  two  sinusoidal  amplitude  distortions  of  0.5 
dB  each,  the  resultant  distortion  needs  to  be  divided  by  2  in  order 
for  it  to  be  0.5  dB . 

From  the  test  run  we  get  both  the  first  and  second  echo 
peak  to  be  -36.2  dB .  The  calculated  first  and  second  echo  peak  is: 


20  lQ9io  + 


10  [%£] 


=  0.0148 


20  log  10  a 


=  36.58  dB 


Tests  4  through  6  are  the  same  as  Tests  1  through  3  respectively, 
except  the  peak  amplitude  distortion  is  set  at  1  dB  instead  of 
0.5  dB  (then  K  =  2).  Test  7  is  the  same  as  Tests  4  or  5  except 
the  distortion  is  introduced  into  the  eighth  echo  pair.  The 
results  are  given  in  Table  4-1. 

TEST  8:  Distortion  is  introduced  in  the  first  and  eighth  echo  pairs 
with  an  amplitude  ripple  of  1  dB  (K  =  2) 


I  =  (10 


sin  ( AX, N ) -1  sin  (AX,N)-1 

*  xin  _ 


where  AY  *  •  AX  =  ^0 

wnere  a  1  215»  a*2  215 


and  0  <  N  <  21 4 


CASE  II;  A  sinusoidal  phase  distortion  of  various  peak-to-peak 
values.  As  defined  previously  the  gain  and  phase  are: 

! 6 |  =  10  log]0  /  I2  +  Q2 
9  -  tan " ^  ^ 


Now  we  want  to  have  a  sinusoidal  phase  variation,  so  0  *  A  sin  (AXN)k. 
To  simplify,  let  us  set  |G|  *  0  dB  so  that: 
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TABLE  4-1.  TEST  RESULTS  OBTAINED  WITH  DEC-20  COMPUTER 


Test  1 
Test  2 
Test  3 
Test  3 
Test  4 
Test  5 
Test  6 
Test  6 
Test  7 
Test  8 
Test  8 


♦Actual  Distortions  outputted  by  the  algorithm  are  in  echo  pairs. 
The  value  shown  is  an  average. 
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1  /  0  \  /  ? 

(1)  10  log10  [I2+Q2]  =  0  and  [I2+Q2]  =  1  or  I2+Q2  =  1 

(2)  j  =  tan  (A  sin(AXN)-k) 

Q  =  I  tan  (A  sin  (AXN)-k) 

Let  A  sin  (  XN)-k  =  0  so  that  Q  =  I  tan  0 
From  ( 1 )  we  have : 

I2  =  1  -Q2 

From  (2)  Q2  =  I2  tan2  0  =  (1-Q2)  tan2  0  =  tan2  0  -Q2  tan2  0 

p  p  9  p 

Q  +  Q  tan  0  =  tan^  0 

Q2  (1+tan2  0)  =  tan2  0 

q  =  tan  0 

1+  tan2  0 

?  1/2 

1  *  rl-  tan^  0i  _  1 _ 

1+tan^  0  (1+tan  9)  c 

0  *  tan”1  ^  =  phase 

To  simplify,  let  us  always  work  in  the  I-V  quadrant  so  that  Q  <  0 

- » — 1 

and  I>^0,  therefore  when  tan  0  <_  0,  then  /  1  +  tan  0  >_  0.  In  order 
to  fit  the  RADC  format  ( 2 ' s  complement)  we  multiply  I  and  Q  by  128. 

TEST  1:  Sinusoidal  phase  distortion  of  6°  peak  value  is  intro¬ 
duced  in  the  first  echo  pair 

0  *  3  sin  ( AXN )  -k  where  AX  =  0  <  N  <  214  and 

k  *  offset  *  42°  places  0  in  fourth  quadrant 


4-3C 


and  I 


q  =  tan(3  sin  ( a x N ) - 4 2 ) 
^1+  tan2(3  sin(AxN)-42) 


/l  -o2 


From  computer  simulation  we  get  the  first  echo  to  be  -31.49  dB . 
Analytical ly  we  get 

2 tt x  b, 

20  1og10  (4lT360-)  =  ech0  in  dB 
where  b^  =  ripple  in  degrees 

b,  -  6°,  20  1°9,0  (jiffj)*  -31.64  dB 

TEST  2:  Cosinusoidal  phase  distortion  of  6°  peak  value  is  intro¬ 
duced  in  the  first  echo 


0  =  3  cos  ( AXN )  -42 


The  computer  simulation  is  basically  the  same  except  that  for  an 
odd  (sin)  distortion  function  the  phase  angle  for  the  two  echoes 
is  180°  apart. whereas  for  an  even  (cos)  distortion  function  the 
two  echoes  are  in  phase.  The  opposite  is  true  for  amplitude  dis¬ 
tortion.  That  is,  for  odd  function  distortions  the  phase  angle 
for  the  two  echoes  is  in  phase  whereas  for  even  distortion  function 
the  phase  angle  for  the  two  echoes  is  180°  apart. 

TEST  3:  Sinusoidal  phase  distortion  of  22.92°  peak  value  is 
introduced  in  the  first  echo 

0  =  11.46  sin  (AXN)  -42 

Calculated  result  is  20  log.|Q  ( ^%x'36Q~‘'")  =  _2B  dB  and  the 
simulated  result  is  -20.02  dB. 


The  results  for  the  phase  distortion  test  cases  are 
shown  in  Table  4-2. 
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TABLE  4-2.  TEST  RESULTS  OBTAINED  WITH  DEC-20  COMPUTER 


Peak  Phase 

Di stortion 
( deqrees  ) 

Echo 

Simulated 

Di stortion 
Phase  (dB) 

Cal cul ated 
Di stortion 
Phase  (dB) 

Test 

1 

6 

1 

-31  .49 

-31  .64 

Test 

2 

6 

1 

-31  .49 

-31 .64 

Test 

3 

22.92 

1 

-20.02 

-20.0 

4.5  PROGRAM  TEST  ON  THE  HP  2100A 

Three  tests  were  run  on  the  algorithm  contained  in  the 
HP  2100A  computer  at  RADC.  The  first  test  is  a  single  echo  test 
involving  a  peak  sinusoidal  amplitude  distortion  of  0.5  dB.  The 
second  test  is  a  double  echo  test  involving  a  peak  sinusoidal  amp¬ 
litude  distortion  of  0.5  dB  with  one  sine  wave  at  twice  the  fre¬ 
quency  of  the  other.  The  third  test  is  a  single  echo  test  involv¬ 
ing  a  peak  phase  distortion  of  six  degrees. 

To  perform  the  first  test,  a  program  has  been  implemented 
(now  on  disk  cartridge  45)  called  "AT10".  A  listing  of  this  program 
...is  gi^en  in  Appendix  3. 

AT10  outputs  M  sample  points  of  a  distorted  sine  wave 


given  by: 


I  =  10 


[sin  AXN]-1 
40 


where:  AX  =  and  N  =  0,lf2,...M 

In  our  case,  we  have  run  AT10  with  256  sample  points  (M  =  255). 
As  mentioned,  the  phase  comoonent  of  this  test  is  0.  AT10  pro¬ 
duces  both  the  single  echo  pair  and  two  pair  test  cases.  The 
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case  just  described  (one  echo  pair)  is  printed  out  in  a  data  file 
called  "AD000";  this  file  contains  256  octal  mode  (06  integers, 
in  two  records,  128  numbers  each;  this  file  is  the  sequential 
output  of  the  "I"  or  amplitude  function  as  N  goes  from  1  to  256. 

To  run  the  double  echo-pair,  amplitude-distortion-only 
case,  we  executed  ATI  0  again  but  such  that  AT10  adds  to  the  first 
waveform,  a  second  one  given  by 

sin(AX2N)-i 


where:  AX2  =  --y  and  N  =  0,1,2 . M  =  255 

The  superimposed  waveform  (this  is  what  is  printed  out  into  the 
data  file)  is 

'sup  ■  O  *  >2>/2 

Again,  there  are  256  samples,  written  into  two  128-number  records 
in  06  format.  The  file  containing  ISup  is  called  "AD001".  This 
double  echo-pair  that  is  equivalent  to  Test  3  of  Table  4-1  which  shows 
the  theoretical  peak  distortion  to  be  -36.58  dB . 

For  the  third  case,  phase  distortion  only,  the  program 
"TEST  1"  was  executed.  The  waveform  produced  by  "TEST  1"  is  a  sine 
whose  phase  is  given  by: 

tan  (3  s in(  AX]  N  )-42) 

q  =  -  >  x  .  360  N  =  oj  >2,. ...255 

1+tan*  (3  sin  (AXm)-42) 

To  run  this  test  case,  we  executed  "TEST  1"  (also  on 
disk  45),  and  used  the  data  file  "AD002".  This  results  in  a  test 
case  of  6°  phase  distortion  equivalent  to  Test  1  in  Table  4-2. 
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Tables  4-3  through  4-5  give  printouts  of  the  data 
contained  in  the  output  files  for  these  three  test  cases,  "AD000", 
"AD001",  and  "AD002"  respectively. 

We  can  now  compare  the  results  obtained  with  the  test 
cases  run  on  the  DEC-20  and  HP  2100A  computers  with  their  respec¬ 
tive  equalizer  algorithms.  Slight  differences  are  to  be  expected 
because  of  the  differences  in  precision  and  the  internal  mathematical 
algorithms.  Table  4-6  shows  the  results. 

The  comparison  shows  excellent  agreement  between  the 
two  equalizer  algorithm  outputs  (DEC-20  and  HP  2100A)  and  demon¬ 
strates  that  they  are  performing  the  way  they  are  supposed  to. 
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TABLE  4-3.  OUTPUT  TEST  1  (ONE  ECHO  PAIR) 
File  ADOOO 


Tap  No. 

Attenuator 

(dB) 

(Turns) 

Line  St 
(deg) 

retcher 

(Div) 

8A 

- .  59754E+02 

9.924 

.  77894E  +  03 

7.2 

7A 

- .  57  500E+02 

9.133 

.  84826E  +  03 

7.8 

6A 

- .  50899E  +  02 

9.931 

.72472E+03 

6.7 

5A 

- . 54626E+02 

9.280 

.1 2426E+04 

1 1  .0 

4A 

- .  7677  5E+02 

9.870 

.1391 6E+03 

1  .3 

3A 

- .  57659E+02 

9.054 

.  48888E  +  03 

4.5 

2A 

- .  63654E+02 

8.21  3 

.  82400E  +  03 

7.6 

1  A 

- .  30575E+02 

2.894 

.  94832E  +  03 

8.8 

0 

.  95424E+01 

0.000 

.  44309E  +  03 

4.1 

IB 

- . 301 02E+02 

2.571 

.4051 9E+02 

0.4 

2B 

- .  581  78E+02 

8.552 

.  60589E+03 

5.6 

3B 

- . 53561 E+02 

9.354 

.  47294E+03 

4.4 

4B 

- . 57641 E+02 

9.961 

. 801 04E+03 

7.4 

5B 

- . 52656E+02 

9.388 

. 1 5492E+03 

1  .4 

6B 

- . 68077E+02 

9.716 

. 23711 E+03 

2.2 

7B 

-.55111 E+02 

9.231 

. 61 405E+03 

5.7 

8B 

- . 58767E+02 

9.938 

.91 032E+03 

8.4 
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TABLE  4-4.  OUTPUT  TEST  2  (DOUBLE  ECHO  PAIR) 

File  ADOOl 


Tap  No. 

Attenuator 
( dB) 

(Turns) 

Line  Stretcher 
(deg) 

(Div) 

8A 

- . 58976E+02 

9.646 

.78231 E+03 

7.2 

7A 

- . 58854E+02 

9.603 

. 86500E+03 

8.0 

6A 

- .  50607E  +  02 

9.816 

.8591 9E+03 

7.9 

5A 

- .  55226E+02 

9.481 

. 90270E+03 

8.4 

4A 

- .  74326E  +  02 

9.505 

. 32962E+03 

3.0 

3A 

-  .  56551  E  +  02 

8.710 

.361 37E+03 

3.3 

2A 

-  .36679E+02 

3.518 

. 93857E+02 

8.7 

1A 

- . 36044E+02 

3.898 

. 9501 2E+03 

8.8 

0 

. 95424E+01 

0.000 

.  44309E+03 

4.1 

IB 

- . 35637E+02 

3.589 

.  39872E+03 

3.7 

2B 

-. 35844E+02 

2.950 

.  491  32E+03 

4.5 

3B 

- . 5281 0E+02 

9.094 

.  60045E+03 

5.6 

4B 

-  . 5731  6E+02 

9.856 

.  97058E+03 

9.0 

5B 

- . 531 65E+02 

9.563 

.  1 3478E+03 

1  .2 

6B 

- . 66063E+02 

9.314 

.  1 0264E+03 

0.9 

7B 

- .  561 90E+02 

9.654 

.59731 E+03 

5.5 

8B 

- .  58025E+02 

9.706 

.  90695E  +  03 

8.4 
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TABLE  4-5.  OUTPUT  CASE  3  (PHASE  DISTORTION  ONLY) 

File  AD004 


Tap  No. 

Attenuator 

(dB) 

(Turns ) 

Line 
( deg) 

Stretcher 

(Div) 

8A 

- . 59834E+02 

9.953 

.  79680E+03 

7.4 

7A 

-.59511E+02 

9.837 

.  73877E+03 

6.8 

6A 

-  .  50643E  +  02 

9.831 

.  96246E  +  03 

8.9 

5A 

- . 56282E+02 

9.841 

.11 579E+04 

10.7 

4A 

- . 72772E+02 

9.269 

.  38686E+03 

3.6 

3A 

- .  57695E  +  02 

9.066 

.  39645E  +  03 

3.7 

2A 

- .  66269E+02 

8.720 

.  67455E+03 

6.2 

1A 

-.311 02E+02 

2.984 

.  1 0384E  +  04 

9.6 

0 

.  95424E  +  01 

0.000 

.44309E+03 

4.1 

IB 

-.31 104E+02 

2.742, 

. 1 7235E+03 

1  .6 

2B 

-  .  58022E  +  02 

8.500 

.  44508E  +  03 

4.1 

3B 

- .  53881 E+Q2 

9.466 

. 78023E+03 

7.2 

4B 

- .  57201  E+02 

9.818 

.8651 8E+03 

8.0 

5B 

- .  53831  E+02 

9.794 

.  1  0693E+03 

1  .0 

6B 

- .  6851  3E+02 

9.803 

. 20575E+03 

1  .9 

7B 

- .  56951  E+02 

9.961 

.  33624E+03 

3.1 

8B 

- . 57922E+02 

9.674 

.  66682E+03 

6.2 
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TABLE 

4-6. 

COMPARISON  OF  TEST  RESULTS  FOR  THE 
EQUALIZER  ALGORITHMS 

DEC-20  AND 

HP  2100A 

Test 

Type  of 

Di  stortion 

Echo 

Affected 

S imul ated 
Level  in 
HP  2100A 

Echo 

dB 

DEC-20 

Theoretical 
Level  in  dB 

Test 

1  , 

0.5  dB 

1  A 

-30.58 

-30.44 

-30.57 

One 

Echo 

Pair 

Peak 

Amp! i tude 

Onl  y 

IB 

-30.10 

(AVG) 

-30.57 

Test 

2, 

1A 

-36.04 

36.36 

-36.58 

Doubl e 
Echo 

Pair 

0.5  dB 

IB 

2A 

-35.64 

-36.68 

(AVG) 

-36.58 

-36.58 

2B 

-35.84 

-36.58 

Test 

3, 

6°  Peak 

1A 

-31  .10 

-31 .49 

-31 .64 

Phase 

Distortion 

Only 

IB 

-31  .10 

(AVG) 

-31 .64 
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While  this  is  not  100%  conclusive,  it  does  provide 
a  high  degree  of  confidence  that  the  algorithm  is  working  properly 
The  only  real  test  is  the  use  of  the  algorithm  over  many  real 
situations  which  completely  exercise  it.  Some  comments  about 
what  to  expect  are  in  order. 

First  of  all,  the  Microwave  Transversal  Equalizer 
(MTE)  to  be  used  to  provide  corrections  to  the  distortion  has 
its  own  iherent  second  order  distortions  which  are  not  constant. 
That  is,  under  one  set  of  amplitude  and  phase  adjustments  for 
each  of  the  taps  there  is  a  given  "self-distortion"  which  is 
included  in  the  data  supplied  to  the  algorithm  to  analyze. 

When  the  amplitude  and  phase  settings  are  now  changed  to  cor¬ 
rect  for  the  total  distortion  reflected  in  the  data,  then  the 
new  settings  of  the  MTE  will  correct  much  of  the  distortion  but 
because  the  MTE  settings  are  different  from  the  original  settings 
there  will  be  a  new  "sel f-di stort i on ”  introduced  by  the  MTE. 

This  means  that  the  correction  is  imperfect  and  the  process  must 
be  iterated  until  the  distortion  residue  is  less  than  the  allow¬ 
able  value  (or  until  the  "self-distortion"  changes  are  greater 
than  the  corrected  indicated  by  the  algorithm). 

Another  source  of  distortion  lies  in  the  hardware 
used  to  provide  the  data  samples.  These  distortions  are: 

1.  the  phase  response  of  the  components 

2.  the  amplitude  response  of  the  components 

3.  noise  introduced  on  the  data  (Gaussian) 

4.  random  spikes  (transients  induced  by 
i nterference ) 
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5.  sinusoidal  variations  (coupled  interference 
from  motors,  oscillators,  etc.) 

6.  timing  errors  (clocks,  circuit  switching 
speeds,  etc.) 

7.  inaccuracies  in  data  samples  (improper 
seal i ng ,  etc  . ) 

8.  digitizing  errors  (bit  errors  in  digital 
operations) 

These  errors  introduce  their  own  "self-distortion"  and  just  as  in 
the  MTE  limit  the  residual  distortion  that  can  be  obtained.  In 
addition,  some  of  the  errors  are  random  so  that  iterations  won't 
help. 

Finally,  there  is  the  algorithm  itself.  It  cannot  be 
considered  fully  debugged  until  all  combinations  of  operation  have 
been  exercised.  This  is  not  a  practical  thing  to  do  in  the  labor¬ 
atory  but  instead  requires  a  perigd  of  field  use.  In  addition,  the 
algorithm  has  limits  of  precision  which  will  ultimately  limit  the 
amount  of  distortion  that  can  be  corrected  even  if  the  hardware 
were  perfect,  although  this  is  the  area  that  is  least  likely  to 
present  a  problem. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


Computer  algorithms  have  been  successfully  developed  to 
provide  open-loop  adaptive  control  of  the  MTE .  The  algorithms  are 
based  upon  the  application  of  FFT  techniques,  and  the  necessary 
corallary  software  programming  procedures  have  been  developed  and 
described . 

Verification  of  the  aforement i oned  FFT  and  associated 
software  programming  has  also  been  accomplished  at  AIL  with  a 
DEC-20  computer.  Such  verification  has  involved  the  aeneration 
of  output  adjustment  data  for  MTE  control  based  upon  computer 
analysis  of  artificially  simulated  time  sidelobe  input  distortion 
levels.  A  similar  procedure  will  be  undertaken  at  RADC  using 
real-time  I  and  Q  data  acquired  at  the  A.  Froelich  High  Power 
Tube  Facility.  Such  distortion  levels  will  be  processed  at  RADC 
using  translation  of  format  programs  between  the  AIL  DEC-20  and 
the  RADC  HP  2100A  computers.  Preliminary  data  suggests  that  the 
post-delivery  program  translation  will  be  successful  after  the 
normal  debugging  procedures  have  been  completed. 

AIL  recommendations  for  future  MTE  related  efforts 
include  the  following: 

«  Continuation  and/or  extension  of  contracted 
efforts  to  include  AIL  post-delivery  support 
for  RADC. 
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•  Investigation  and  definition  of  suitable  closed- 
loop  techniques  and  procedures  for  adaptive  oper¬ 
ation  of  MTE . 

•  Development  of  suitable  solid-state  time  delay 
devices  to  replace  the  present  manually  adjust¬ 
able  line  stretchers. 

•  Modification  of  the  MTE  to  expand  capability  by 
providing  32  tap  operation. 

•  Development  of  suitable  electronic  interface  to 
implement  fully  adaptive  closed-loop  operation 
of  the  MTE. 
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APPENDIX*  1 

SAMPLE  COMPUTER  RUN  WITH 
INPUT  DATA 
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APPENDIX  1:  SAMPLE  RUN 


This  appendix  includes  a  sample  run  of  the  program  with  the  input  data 

entitled  AD000. 


:  PU»  RIL 
1 

INITIRL  SETTINGS  RT  STEP  0  SHOULD  BE: 


RTTENURTQP 


LlH.E  STRETCHER 


1  R= 10. 000 

2  H=1 0. 000 

3  R=1  0.000 

4  R*1  0.000 

5  ft=l  0.000 

6  R=1 0.000 
?  R*1  0.000 
3  R=1 0.000 

M  *  0.000 

1  B=1  0.000 

2  B*1 0.000 

3  B=1  0. 000 

4  B=1  0.000 

5  B*1 0.000 

6  B=1 0.000 

7  B=1  0.000 
3  B*1  0.000 


l 

R= 

7.9 

c! 

R  = 

6. 2 

3 

R= 

*  * 

4 

H  = 

1.5 

cr 

..i 

R= 

9.  1 

Initial 

6 

R= 

7.5 

Attenuator 

r’ 

H= 

7.  0 

Setti nos 

3 

R= 

6.  1 

M  = 

4.  1 

1 

B= 

1.2 

c! 

B= 

3.7 

3 

B= 

5.2 

4 

B= 

7.2 

5 

B= 

.5 

6 

B= 

1.4 

7 

B- 

3.2 

3 

B= 

6. 2 

Initial 

Line 

Strecher 

Settings 
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APPENDIX  1  SAMPLE  RUN  (continued 
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APPENDIX  1  SAMPLE  RUN  (continued) 


INTERMEBIRTE  SETTINGS  RT  STEP  1  BRE: 

RTTENURTOR  LINE  STRETCHER 


1 

R= 

1 . 976 

1 

R= 

O  ^  C; 

a 

R= 

6.  a  03 

3 

fi= 

i  ■  Ca 

3 

R= 

9.319 

3 

R= 

8.  9 

4 

R= 

8 . 486 

4 

R= 

1 . 6 

5 

R* 

9 . 669 

5 

R= 

8. 6 

6 

R= 

9.705 

6 

R= 

3.  9 

Ra 

9.  338 

7 

R= 

6  -  3 

3 

R= 

9.434 

3 

R= 

6*  1 

M  = 

0.  000 

n  = 

4.  1 

1 

£= 

1.573 

i 

£= 

.4 

c 

B- 

6 .  085 

C. 

B= 

5-5 

3 

B» 

9.  665 

3 

B* 

6.  0 

4 

B~ 

9.  363 

4 

B= 

: .  i 

5 

B* 

9.  707 

5 

B= 

A  •  0 

6 

B= 

3.  934 

6 

B= 

w*  ■ 

B= 

9.  343 

“7 

B* 

3.  9 

s 

B= 

9.503 

3 

B= 

8  •  c! 

RFTER  0  I TERRTIONS  DISTORTION  CRN  NDT,£E  COMPENSRTED 
RIL  :  STOP 


Note:  XD  *  20  in  $AT10  when  generating  input  for  this  run. 

In  subsequent  runs,  XD  *  40 
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APPENDIX  2 

PROGRAM  TO  GENERATE 
TEST  DATA 
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APPENDIX  2.  PROGRAM  TO  GENERATE  TEST  DATA 
This  routine  generates  ideal  sinusoidal  distortion  as  the  data  file  for  the  Dreary 


0  0  0 1 

FTN4. 

L 

00  Oil' 

PROGRAM  hTIO 

0  0  0  3 

DIMENSION  I BUF 8 < 358 > *  I BUF 3 < 858' « MBUF <858> 

0  0  04 

D I  MEN  I  I  ON  1  DC  £  •  878  •  .  N*M  <  3 •  I  BUF  0  <85 81'  .  I  SIZE  <8 

0  0  05 

D I  MEM  I  I  ON  LBUF  •  858  -'  .  I  BUF  1  <  858  ' 

00  0*3 

D I  MEN  3  ION  I  DC  S3  <878  •  •  NAM1  <..3> 

i  j  0  0  7 

I  TYPE =8 

0  0  08 

I3IZE=8 

0  0  03 

I  SIZE <8>  =183 

0  0 1  0 

DhTR  NFtM  ■"  8HhD*  8H00«  8H0 

0  0 1  1 

DhTh  MhM  1  -'8HhD«  8N00*  8H1 

0  0 1  8 

ZhLL  CR'EFtI  IDCB*  IEFR.NRM.  I3IZE*  ITYFE) 

>:•  o  i  ? 

ChLL  CR'EFtT  <  I  DC  B ! «  I ERP • N*M1 • I SIZE* I T  YPE> 

0  0 1  4 

N=858 

0  0 1 5 

X  D  =  4  0 . 

0  0  l  8 

DX1 =380.  • FLOhT  <N-1 • 

0  'J  1  7 

DX8=DX1>8. 

u  0 1  3 

DO  3  0  .1=1*  N 

0  0 1 8 

X=D:<1  ♦FlOhT  •:  J-l 

0  08  0 

Y=D7.8  *Fi_OhT  0-1  1 

0  0  8  1 

l  =  *  i  o .  ■  <:i m  < :<:♦  3 .  i4i 5*  1 3 o .  > - 1  > /xio  > 

0  0.38 

X8=<  1  0.  ♦♦<  <' S I N  <  Y ♦  3 .  1415  *180.  >  ^XD>  > 

0  08  3 

LBUF  ■  J'=183»X1 

0084 

XM=.XH-X8‘  8. 

0085 

MBUF > J>=183^XM 

0  088 

I F  ■ MBUF <  J  > . GE . 1 83 ■  MBUF <J> *  187 

0  087 

MBUF  i  J".  =858>MEUF  J> 

0083 

IF  <LBUF  <  J  •  .  EQ.  18S>  LBUF  <1.0  =187 

0083 

LBUF  <  J<  =858»LBUF  <. J> 

0  0  3  0 

80 

CONTINUE 

0  0  3 1 

DO  30  1=1*188 

0  0  38 

i  buf  o  a  •  =lbuf  a  - 

0  0  3  3 

I BUF 8 < I ' =MBUF <  I  > 

00  34 

30 

CONTINUE 

00  35 

DO  4  0  k  -  1  *  1  88 

0  0  38 

I  BUF  1  <  h:  >  =LBUF  <  1  83+K  .• 

0  0  3  7 

I  BUF  3  *  k  =M  BUF  <.  1  83+k.  > 

0  033 

4  0 

CONTINUE 

00  39 

C  hLL  MR  I  TF  <  I  DC  B  *  I  ERR «  I  BUF  0  • 

0040 

C  hLL  MR  I  TF  <  1  DCE  *  I  ERR  -  I  BUF  1  :■ 

004  1 

C HLL  MR  I TF  < I  DC  El* I ERR* I BUF 8 > 

0048 

C  FILL  MR  I  TF  <  I  DC B 1 «  I  ERR  •  I  BUF  3 > 

004  3 

c fill  clo:e<idcb*  ierr- 

0  044 

CFiLL  CLOSE  >IDCB1.IERR' 

0045 

STOP 

0048 

END 
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This  program  generates  a  sinusoid  with  phase  distortion  only. 


LIST?  $TEST 1 : : 45 

tTESTl  T  =  0 0 0 04  IS  ON  CR00045  USING  00005  ELKS  R=0034 
RVL3L1  FTN4 


000  £ 
0  0  0  3 
0  0  04 
0005 
0  0  06 


PROGRAM  TEST 1 

INTEGER  I DCE C£7£)  ?NAMC3) ?  I  BUR  0  (£56) ? IEUF1 (£56) ? LEUF <£56) 
INTEGER  I SIZE <£> 

I TYPE=£ 

N=£56 


0007 
0008 
00  09  C 
■'  0 1  0 
0  0 1 1 
0  0 1  £ 

0  0 1  3 
0  0 1 4 
0  0 1  5 
0  0  i  6 
0  0 1  7 
0  0 1 3 
0  0 1  9 
0  0£  0 
00£1 
0  0££ 

0023 
00£4 
00  £5 
0  0  £  8 
0  0£7 
00£3 
00  £9 
0  0  3  0 
0  0  3 1 
0  0  3£ 

0 1 1 3  3  1  0 

00  34 
8035 
0038 

0  0  37  3 1 } 

00  38 
0039 
004  0 
0041 


ISIZE=£ 

ISIZE  C£> =1£8 

CALL  CREPT  Cl  DCE?  I  ERR?  NAM?  ISIZE?  I  TYPE) 
DATA  NAM/  £HAD? £H00? £H4  / 

CFiLL  CREPT  Cl  DCE?  IERR?NAM?  ISIZE?  I  TYPE) 
DX 1 =36 0 .  FLOAT CN- 1 > 

DO  10  1=1 ?N 

xi=dxi ♦float  <  i-n 

X=  3.  ♦  S I N  (X 1  ♦  £ .141 5/  ISO.)  -4£. 

X£=S  I N  <»3 .  1415--18  0 . ) 

X3  =  CDS  <X»3. 1415/180. ) 

T=X£/X3 
T £=T  ♦  T 

X4=SQRT  Cl+T£> 

IFCT.GT.O)  X4=-X4 
Q=T/X4 
I Q= 1 £8 . ♦Q 
XI=SQRT<1.-Q^Q) 

1 1  =  1 £8^X1 

IFCII.GE. 1£7)  1 1  =  1 £7 
IF  <  II.LE.-1£8>  1 1 =-l£3 
IF  (IQ. GE. 1£7>  I Q= 1 £7 
IF  CIQ.LE.-128)  IQ=-1£8 
IFCII.LT. 0)  II=II+£56 
IF  CIQ.LT.O-  IQ=IG+£56 
LBUF  Cl) =£56»I I  +  IQ 
CONTINUE 
DD  30  J= 1 ? 1 £8 
IB UFO CJ) =LBUF CJ) 

I EUF 1 (J> =LEUF CJ+128) 

CONTINUE 

CALL  WR I TF  C I DCE ? I ERR ? I EUF 0) 

CALL  WRITFcIDCB? I ERR? IEUF1) 

CRLL  CLOSE  Cl DCE? I ERR) 

STOP 
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Input  for  Case  1:  Octal  Format 


:LIST» ADOOO 

RIi OOO  T=  0000c-  IS  ON  CROOOO 


REC«  00001 

074000  074400  074400  074400 
075000  075000  075000  075000 
075400  075400  075400  076000 
076000  076000  076400  076400 
076400  076400  077000  077000 
077000  077000  077000  077400 
077400  077400  077400  077400 
077400  077400  077400  077400 
077400  077400  077400  077400 
077400  077400  077400  077400 
077400  077400  077400  077400 
077000  077000  077000  077000 
076400  076400  076400  076400 
076000  076000  076000  076000 
075400  075400  075400  075000 
075000  074400  074400  074400 

RECjs  000  02 

074000  074000  074000  074000 
073400  073400  073400  073000 
073000  072400  072400  072400 
072000  072000  072000  072000 
071400  071400  071400  071400 
0714  0 0  0714  0  0  0  710  0 0  0  710  0  0 
07100  0  0  7 1 0  0  0  0  710  0  0  0  7100 0 
0710 0  0  071 0 1 j  0  0  7 1 0  0  0  0  7 1 0  0  0 
071000  0710 0 0  07 1 0 0 0  07 1 0 0 0 
071000  0710  0  0  0  7 1 0  0  0  0  7 1 0  0  0 
071000  0  7 1 0  0  0  0  7 1 0  0  0  0710 0  0 
071400  071400  071400  071400 
072000  072000  072000  072000 
072400  072400  072400  072400 
073000  073000  073000  073000 
073400  073400  073400  074000 


USING  00002  BLkS  R=0128 


074400  074400  074400  075000* 
075000  075400  075400  075400* 
0 760  0  0  0 7 6  0  0  0  0 7 6  0  0  0  076 0  0  0  * 
076400  076400  076400  076400* 
0  7  7  0  0  0  0  7  7  0  0  0  0  7  7  0  0  0  0  7  7  0  0  0  ♦ 
077400  077400  077400  077400* 
077400  077400  077400  077400* 
077400  077400  077400  077400* 
077400  077400  077400  077400* 
077400  077400  077400  077400* 
077400  077400  077000  077000* 
077000  077000  077000  076400* 
076400  076400  076000  076000* 
076000  075400  075400  075400* 
075000  075000  075000  075000* 
074400  074400  074400  074000* 


074000  073400  073400  073400* 
073000  073000  073000  073000* 
072400  072400  072400  072400* 
072000  072000  072000  072000* 
071400  071400  071400  071400* 
0 7 1  0  0  0  0 7 1  0  0  0  0  7 1  0  0  0  0  7 1  0  0  0  ♦ 
0  7 1  0  0  0  0  7 1  0  0  0  0 7 1  0  0  0  0  7 1  0  0  0  ♦ 
0  7 1 0  0  0  0  7 1 0  0  0  071 0  0  0  0  7 1 0  0  0 ♦ 
071 0  0  0  0  710  0  0  0  710  0  0  071 0  0  0 ♦ 
0  7 1 0  0  0  0  7 1 0  0  0  071 0  0  0  071 0 0  0 ♦ 
0  710  0  0  0714  0  0  0714  0 0  0714 0  0 ♦ 
071400  071400  071400  072000* 
072000  072000  072000  072400* 
072400  072400  073000  073000* 
073000  073400  073400  073400* 
074000  074000  074000  074000* 
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Input  for  Case  2:  Octal  Format 

LIST*  FiD  001 

RDO 0 1  1=000 OS  IS  ON  CP00002  USING  0000c:  BLKS  P=0128 

PEC"  00001 


074000 

074400 

074400 

074400 

075000 

075400 

075400 

075400 

076000 

0  7 6  0 :  j  0 

076400 

076400 

077000 

077  0 0 0 

077  000 

077  0  0  0 

077000 

077  0 0 0 

077400 

077400 

077400 

0774 0 0 

077400 

077000 

077000 

077000 

077000 

077000 

076400 

0764  0 0 

076400 

076400 

076000 

076000 

076 0  00 

075400 

075  0 0  0 

075 0  0 0 

C 7  5  0  0  0 

075000 

074400 

074400 

074400 

074400 

074000 

074000 

074000 

074000 

073400 

073400 

073400 

073400 

073400 

073400 

073400 

073400 

073400 

073400 

073400 

073400 

074  0 0 0 

074000 

074000 

074  000 

PEC-‘:  0000£ 

074000 

074000 

074400 

074400 

074400 

074400 

074400 

074400 

075000 

075  0  0  0 

075000 

075000 

075  0 0  0 

075000 

075000 

075000 

074400 

074400 

074400 

074400 

074400 

074400 

074 000 

074000 

074000 

073400 

0734  0 0 

073400 

073000 

073000 

073000 

073  0  0  0 

07 £4 00 

072400 

072400 

072000 

072000 

0714  0  0 

071400 

071400 

071400 

0714  0  0 

071400 

071400 

071 000 

071 000 

071 000 

071 000 

071400 

0714  0  0 

071400 

071400 

071400 

072000 

072000 

072000 

072400 

072400 

0724 0 0 

072400 

073400 

073400 

073400 

073400 

0744 on 

075000 

075000 

075000* 

0754 00 

076000 

076000 

076000* 

076400 

0764  0  0 

076400 

076400* 

077000 

077000 

077000 

077000* 

077400 

077400 

077400 

077400* 

077000 

077000 

077000 

077000* 

077000 

077000 

0764 00 

076400* 

0764  0  0 

076000 

076000 

076000* 

075400 

075400 

0754 0 0 

075400* 

075000 

075000 

074400 

0744  0 !:* 

0744 0  0 

074000 

074000 

074  000* 

074000 

074000 

074000 

074 OuO* 

073400 

073400 

073400 

073400* 

073400 

073400 

073400 

073400* 

073400 

074000 

074000 

074000* 

074000 

074000 

074000 

074000* 

074400 

074400 

074400 

0744  0  0* 

0744 0  0 

074400 

074400 

075000* 

075000 

075000 

075000 

075000* 

075000 

075000 

075000 

075000* 

074400 

074400 

074400 

074400* 

074000 

074000 

074000 

074000* 

0734 0  0 

073400 

073400 

073000* 

073 0  0  0 

072400 

072400 

0724 0  0* 

0  72  0  0  0 

072000 

072000 

072000* 

071400 

071400 

071400 

071400* 

071 000 

071 000 

071 000 

071000* 

071 000 

071 0  0  0 

071 000 

071400* 

071400 

071400 

071400 

071400* 

072000 

072000 

072000 

072400* 

073000 

073000 

073000 

073000* 

074000 

074000 

074000 

074000* 
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Input  for  Case  3:  Octal  Format 


LIST*  hD  0 04 

RD004  T=00002  IS  ON  CR00002  USING  00002  ELKS  R=0128 


REC"  00001 

057653 

057653 

057653 

057653 

057654 

060254 

06  0254 

060254 

060255 

06 0255 

06  ft  6  55 

06  i'ik55 

06  0656 

060656 

060656 

060656 

061256 

061257 

061257 

061257 

061 257 

061257 

061257 

061257 

061660 

061660 

06 1 66 0 

06 1 6  60 

06 1 66  0 

061660 

06 166  0 

06 1 66  0 

061660 

06 166  0 

06 166  0 

06 1 66 0 

06 1 66  0 

06 1 66  0 

06 1660 

06 166  0 

06 1 66 0 

061657 

06 1 657 

06 1657 

061257 

06 1 257 

061257 

061257 

06 1256 

06 1  c56* 

061256 

060656 

060656 

060656 

06  0655 

06 0655 

060255 

060255 

06  0254 

060254 

nr  .cme 

057654 

057654 

057654 

057653 

REC»  00002 

057653 

057253 

057253 

057252 

057252 

057252 

056652 

05665 1 

056651 

056651 

05665 1 

05665 1 

056250 

056250 

056250 

056250 

055647 

055647 

055647 

055647 

-BP-  COItE-CPR 

055647 

055647 

055647 

055647 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055646 

055647 

055647 

055647 

055647 

055647 

055647 

056250 

05625  0 

05625  0 

056250 

056250 

056250 

056251 

05665 1 

056651 

05665  1 

05665 1 

057252 

057252 

057252 

057252 

057653 

057653 

057654 

057654* 

06 0254 

060254 

06  0255 

060255* 

06 0655 

060655 

06  0655 

060656* 

060656 

061256 

061256 

061256* 

061257 

061257 

061257 

061257* 

06 1 257 

061657 

061657 

06 1 66 0* 

06 1 66  0 

06 1 6-6-  0 

06 166 0 

u6 1 66 0* 

061660 

06 1 66 0 

06 i 660 

06- 1 66  0* 

06 1 6-6  0 

06 1 66  0 

1J  tlt.6  Li 

06 1 66 0* 

06 1 660 

06 1 66 0 

06 1 66 0 

o  o  1 t*  o  ♦ 

06 1 257 

061257 

061257 

061257* 

061257 

061257 

061257 

1 

060656 

06  0656 

06  0656 

06  0656* 

06  0655 

06  0655 

06  0255 

060255* 

060254 

06  0254 

06  0254 

060254* 

057653 

057653 

057653 

057653* 

057252 

057252 

057252 

057252* 

05665 1 

056651 

056651 

056651* 

056250 

056250 

i.i  5  k  d  5  U 

•55625  0* 

056250 

056250 

056250 

055647* 

055647 

05564  7 

055647 

055647* 

055647 

055246 

055246 

055246* 

055246 

055246 

055246 

055246* 

055246 

055246 

055246 

055246* 

055246 

055246 

055246 

055246* 

055246 

055246 

055246 

055246* 

055647 

055647 

055647 

055647 * 

055647 

055647 

055647 

0556.47  ♦ 

056250 

05625 0 

056250 

05625 0* 

056651 

056651 

056651 

056651* 

056652 

057252 

057252 

057252* 

057252 

057253 

057653 

057653* 
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'--V /;:*  * 


'V 


1 


Input  for  Case  3: 
Decimal  Format 


24491 

24491 

24491 

24491 

24491 

24491 

244?2 

29492 

24492 

24746 

2474S 

24748 

2474?. 

24743 

2474? 

2474? 

24749 

24749 

25005 

2*005 

2505)5 

25005 

2:003 

23006 

2500* 

25006 

15006 

25006 

25006 

1:262 

25262 

25262 

25262  • 

-,c  *.27 

-  t 

-  C  •  2  -3 

25263 

2351? 

2:31? 

:3'7? 

25520 

25520 

7  ”.-70 

25520 

25520 

OCC-fj 

25520 

25520 

25520 

25520 

25*20 

25:20 

-752u 

■3320 

23520 

25520 

25520 

15510 

25520 

25520 

25320 

25520 

25320 

2:320 

2-5 

25520 

25*20 

25520 

23520 

1*520 

25520 

25520 

2551? 

c  r  3 

255  i? 

23263 

15263 

25263 

25263 

25263 

25263 

25263 

3 

25263 

23262 

^  -  7  *’.? 

IC.-j.-ii 

2~  v'6 

";S'V>r 

1!VmV. 

25006 

2500* 

25005 

25005 

23005 

25005 

24749 

24749 

2474? 

2474? 

2474? 

24748 

2474s 

2474,;. 

24748 

24749 

24492 

244?2 

24492 

24491 

244?  1 

24491 

244?1 

24491 

24451 

24235 

24235 

24234 

24234 

24239 

24234 

242:74 

24224 

24234 

23?78 

23?77 

?397? 

23977 

23227 

13-77 

22977 

23?77 

23477 

23  977 

23720 

27720 

23720 

22720 

23720 

23720 

23720 

23720 

23720 

23720 

23720 

23463 

23463 

23463 

274  a.? 

27463 

23463 

23463 

224:3 

234*3 

23463 

23463 

28463 

23463 

23.463 

23206 

•r^p  i 

232u6 

23206 

33206 

27206 

23206 

23206 

23206 

222)6 

4.J4VV 

23205 

23206 

23206 

23206 

23206 

23206 

2:^6 

23206 

23206 

23206 

23206 

23206 

23208 

23206 

wi.VS 

23206 

23206 

23206 

23306 

23206 

23206 

23206 

23206 

23206 

23206 

23462 

23463 

23463 

23463 

23*63 

23463 

23463 

234 6? 

23463 

774.. 3 

23463 

234-3 

28463 

23463 

23453 

23720 

23720 

23720 

23720 

23/20 

22720 

23720 

23720 

23720 

23721 

23977 

2:  ;77 

23977 

23-977 

23977 

23977 

23977 

23?77 

23978 

24234 

24234 

24234 

24234 

24234 

24234 

24234 

24234 

24235 

24491 

24491  STOP 

I 


\ 
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MISSION 
of 

Rome  Air  Development  Center 

MOC  plan*  and  executes  research,  dtveI.opme.nt,  teAt  and 
selected  acquisition  programs  in  support  of  Command,  Control 
Cotmunic.ati.ons  and  Intelligence  (C3I)  activities.  Technical 
and  engineering  support  within  areas  of  technical  competence 
as  provided  to  ESV  Program  Offices  IPOs)  and  other  BSD 
element A.  The  principal  technical  mission  areas  are 
communications ,  electromagnetic  guidance  and  control,  sur¬ 
veillance  of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


—  -r  <pmm mx w >.  i. - 


END 


DATE 
FILMED 


DTIC 


